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Chapter  1 


1.1  Introduction 

An  efficient  closed  form  asymptotic  representation  for  a  grounded  double-layered  anisotropic 
uniaxial  geometry  is  developed.  The  large  parameter  of  this  asymptotic  development  is  di¬ 
rectly  proportional  to  the  lateral  separation  between  the  source  and  observation  point. 
However  the  asymptotic  solution  remains  accurate  even  for  very  small  (a  few  tenths  of 
a  wavelength)  lateral  separation  of  the  source  and  field  points.  The  asymptotic  closed 
form  dyadic  Green’s  function  has  been  cast  in  such  a  form  that  the  physical  behavior  of 
the  electromagnetic  fields  due  to  anisotropy  of  the  medium  reveals  itself  through  a  simple 
mathematical  parameters.  Thus,  the  physical  understanding  of  the  interaction  of  the  spa¬ 
tially  confined  source  with  an  anisotropic  (uniaxial)  double-layered  grounded  slab  is  greatly 
enhanced  through  the  newly  developed  asymptotic  closed  form  representation  of  the  dyadic 
Green’s  function. 

Also,  this  efficient  representation  is  very  useful  in  the  moment  method  (MM)  solution 
of  the  current  excited  on  the  microstrip  antennas  and  arrays  in  a  grounded  double-layered 
uniaxial  geometry,  as  well  as  the  volumetric  current  excited  within  a  dielectric  scatterer 
buried  in  a  grounded  double  layered  anisotropic  uniaxial  slab.  The  MM  analysis,  especially 
for  microstrip  arrays  and  guided  wave  structures,  requires  a  very  large  number  of  compu¬ 
tations  where  the  lateral  distance  between  the  source  and  the  field  points  in  the  dyadic 
Green’s  function  can  range  from  extremely  small  to  very  large  values  [1,2].  The  newly  de¬ 
veloped  asymptotic  Green’s  function  can  drastically  improve  the  computational  speed  and 
the  efficiency  of  the  MM  solution,  as  compared  to  the  use  of  the  conventional  Sommerfeld 
form  of  the  Green’s  function.  For  small  lateral  separation  (less  than  a  few  tenths  of  the 
wavelength)  between  the  source  and  observation  points,  one  can  switch  to  either  modified 
Sommerfeld,  or  to  the  radially  propagating  integral  representation  of  the  dyadic  Green’s 
function  which  also  developed  during  the  project.  For  very  small  lateral  separations,  less 
than  a  quarter  to  half  a  free  space  wavelength,  the  modified  Sommerfeld  representation  is 
more  efficient;  for  the  relatively  small  to  moderately  large  the  radially  propagating  integral 
representation  of  the  dyadic  Green’s  is  should  be  used. 

Hence,  a  systematic  numerical  algorithm  is  developed  which  can  efficiently  be  used  to 
compute  the  dyadic  Green’s  function  for  a  grounded  anisotropic  uniaxial  double-layered 
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slab  for  the  entire  range  of  the  separation  of  source  and  observation  points. 

The  format  of  the  report  is  as  follows:  Chapter  2  deals  with  the  formulation  of  the 
electromagnetic  dyadic  Green’s  function  in  anisotropic  uniaxial  multi-layered  media.  A 
general  formulation  of  the  dyadic  Green’s  function,  and  the  theory  of  2-propagating  modes 
in  anisotropic  uniaxial  media  is  discussed  in  this  chapter.  The  dyadic  Green’s  function 
consists  of  the  integral  over  the  plane- wave  spectrum  of  a  continuum  of  vector- wave  function, 
in  which  the  vector-waves  are  guided  (propagated)  in  the  preferred  z-coordinate  direction. 
Furthermore  a  related  solution  such  as,  the  electromagnetic  dyadic  Green’s  function  for 
an  anisotropic  uniaxial  multi-layered  medium  due  to  periodic  electric  as  well  as  magnetic 
point  currents  has  also  been  developed.  This  solution  is  important  in  many  applications 
involving  anisotropic  uniaxial  stratified  media,  such  as  the  design  of  high  performance  large 
phased  array  in  mulit-layered  envireoment,  the  analysis  of  the  frequency  selective  surfaces 
(FSS),  multi-layered  guided  structures,  and  grating.  In  Chapter  3  a  general  formulation 
presented  in  Chapter  2  will  be  specialized  to  a  canonical  geometry  of  a  grounded  anisotropic 
uniaxial  double-layered  slab  where  the  electric  point  current  located  within  the  slab.  A 
Sommerfeld  and  subsequently  a  radial  propagation  integral  representation  of  the  dyadic 
Green’s  function  of  this  geometry  will  then  be  presented.  Also  in  this  chapter,  two  distinct 
efficient  methods  of  the  numerical  integrations  will  be  presented.  The  first  method  is  based 
on  removal  of  the  singularity  of  the  integrand  of  the  Sommerfeld  (z-propagating)  integral 
representation  of  the  dyadic  Green’s  function.  The  second  method  results  from  the  radially 
propagating  integral  representation  of  the  Green’s  dyadic.  The  uniform  asymptotic  closed 
form  approximation  of  the  dyadic  Green’s  function  for  anisotropic  uniaxial  double-layered 
grounded  slab  is  the  topic  of  Chapter  4.  This  asymptotic  Green’s  function  is  developed 
employing  the  formal  radial  propagation  integral  representation  of  the  Green’s  function 
presented  in  Chapter  3.  Chapter  5  is  devoted  to  numerical  results  based  on  the  accuracy 
and  efficiency  of  our  newly  developed  asymptotic  dyadic  Green’s  function.  The  effect  of 
the  anisotropy  of  the  material  on  the  electromagnetic  field  and  the  physical  interpretation 
of  our  result  will  also  be  presented  in  this  chapter.  The  conclusion  and  the  direction  of 
our  future  research  will  be  presented  in  Chapter  6.  The  analytical  expressions  for  the 
vector-wave  modes,  for  a  general  anisotropic  uniaxial  medium,  is  derived  in  Appendix  A. 
In  Appendix  B  the  reciprocity  and  orthogonality  relations  for  the  vector-wave  modes  are 
presented.  Some  essential  functions  as  well  as  residues  of  the  surface  waves  which  are  being 
used  to  construct  the  radial  propagation  integral  representation  as  well  as  the  uniform 
asymptotic  closed  form  dyadic  Green’s  function  are  presented  in  Appendix  C.  Appendix  D 
devoted  to  the  asymptotic  behavior  of  the  integrands  of  the  Green’s  functions  which  are 
employed  in  singularity  removal  technique. 
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Chapter  2 


Formulation 


2.1  Introduction 

In  this  chapter  a  complete,  plane- wave  spectral,  vector-wave  function  expansion  of  the  elec¬ 
tromagnetic,  electric,  and  magnetic,  dyadic  Green’s  function  for  electric,  as  well  as  magnetic, 
point  currents  for  a  planar,  anisotropic  uniaxial  multi-layered  medium  is  presented.  It  is 
given  in  terms  of  z- propagating,  solenoidal  vector-wave  functions,  where  z  is  normal  to  the 
interfaces,  and  it  is  developed  via  a  utilization  of  the  Lorentz  reciprocity  theorem.  For  the 
completeness  the  electromagnetic  dyadic  Green’s  function  for  periodic  electric  as  well  as 
magnetic  point  current  sources  are  also  presented. 

The  geometry  of  the  problem  is  shown  in  Figure  2.1,  where  the  electrical  parame¬ 
ters,  in  the  layer  are  assumed  to  be  homogeneous  and  uniaxial  (I  =  etIt  + 

ezzz  ;  /x  =  -f  nzzz).  It  is  shown  that  the  electromagnetic  field  at  a  given  point  within 
the  multi-layered  region  consists  of  four  distinct  waves  (two  for  each  vector- wave  type) 
caused  by  the  presence  of  the  multi-layered  medium.  Moreover,  the  presence  of  bound¬ 
aries  in  anisotropic  media  cause  mode-coupling  at  the  planar  interfaces;  these  modes  are, 
in  general,  continuously  coupled  within  the  stratified  region.  However,  due  to  the  reflection 
symmetry  of  the  anisotropic  uniaxial  structure  considered  here  the  modes  are  decoupled, 
and  hence;  the  Green’s  dyadic  within  each  of  the  layers  can  be  expressed  in  terms  of  four 
uncoupled  distinct  waves  (two  for  each  mode),  which  resembles  the  dyadic  Green’s  function 
for  multi-layered  isotropic  medium. 

A  systematic  formulation,  presented  here,  is  useful  in  many  problems  dealing  with  new 
types  of  anisotropic  materials,  such  as  advanced  composites,  ceramics,  and  honeycomb 
structures,  which  are  finding  increasingly  important  applications  in  high-frequency  elec¬ 
tromagnetics  and  optics.  These  applications  range  from  microwave  and  millimeter-wave 
integrated  circuits  and  optical  devices,  to  antenna  radomes  and  radar  absorber  materi¬ 
als.  The  importance  of  this  analysis,  besides  the  physical  significance  of  understanding  the 
interaction  of  electromagnetic  sources  with  complex  media,  is  due  to  a  serious  need  for 
efficient  design  procedures  for  high  performance  monolithic  integrated  circuits  that  operate 
in  the  high  frequency  regime  (millimeter  to  optical  wavelength),  as  well  as  to  a  crucial 
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need  for  an  analytical  model  for  the  non-destructive  evaluation  (NDE)  of  a  variety  of  the 
increasingly-used  advanced  composite  materials. 

Due  to  the  widespread  application  of  complex  materials,  there  has  been  considerable 
research  on  the  interaction  of  electromagnetic  waves  with  multi-layered  anisotropic  media. 
We  restrict  ourselves  to  those  work  which  directly  related  to  anisotropic  uniaxial  media. 
Wait  [3]  and  Kong  [4,  5]  investigated  the  radiation  of  a  dipole  source  over  a  uniaxial  strat¬ 
ified  media.  Interaction  of  the  electromagnetic  source  with  anisotropic  uniaxial  media  has 
also  been  extensively  studied  by  Felsen  and  Marcuvitz  [6].  Ali  and  Mahmoud  [7]  derived 
a  Green’s  dyadic  for  a  uniaxial  stratified  medium,  using  dyadic  Green’s  function  theory 
in  conjunction  with  the  primary  and  secondary  dyadic  Green’s  function.  Tang  [8],  and 
Kwon  and  Wang  [9]  have  investigated  the  problem  of  various  types  of  the  point  sources  in 
uniaxial  stratified  media.  They  derived  the  expressions  for  the  electric  and  magnetic  field 
components  by  satisfying  the  boundary  conditions  at  each  interface  and  a  proper  condition 
at  the  source  point.  Krowne  [10]  has  given  Green’s  function  solutions  for  uniaxial  and  bi¬ 
axial,  layered,  planar  structures  excited  by  a  surface  electric  current.  His  results  are  given 
in  terms  of  three  different  longitudinal-to-an-axis  formulations  of  the  differential  coupled 
field  component  sets.  The  radiation  of  an  electric  dipole  in  the  presence  of  the  anisotropic 
uniaxial  grounded  layered  media  has  been  investigated  by  Tsalamengas  and  Uzunoglu  [11]. 
An  excellent  review  of  wave  propagation  in  general  anisotropic  media  has  recently  been 
given  by  Chew  [12]. 

In  the  present  work  the  complete,  spectral,  vector-wave  expansion  for  the  electromag¬ 
netic  (electric  and  magnetic)  dyadic  Green’s  function  for  electric,  as  well  as  magnetic,  point 
current  dipole  sources  for  planar,  anisotropic  uniaxial  multi-layered  media,  using  solenoidal 
vector-wave  functions  [14]  is  derived.  The  orthogonality  of  the  modes  over  a  plane  surface 
[6]  transverse  to  the  direction  of  the  propagation,  z,  together  with  a  reciprocity  relation  are 
employed  to  construct  the  Green’s  dyadics.  Moreover,  the  method  that  is  employed  here 
utilizes  only  solenoidal  vector-wave  functions  [13];  hence,  the  dyadic  delta  function  term  at 
the  source  point  is  included  explicitly  as  a  correction  to  the  general  solenoidal  vector-wave 
expansion,  which  is  valid  outside  the  source  point  [13,  15]. 

The  format  of  the  chapter  is  as  follows:  Section  2.2  deals  with  the  formulation  of  the 
electromagnetic  dyadic  Green’s  function.  A  general  formulation  of  the  dyadic  Green’s  func¬ 
tion  in  anisotropic  media,  and  the  theory  of  z-propagating  modes  in  anisotropic  media 
will  be  discussed  in  parts  2.2-A  and  2.2-B  respectively.  In  Section  2.3  the  electromag¬ 
netic  dyadic  Green’s  function,  H(r,  r'),  for  electric  and  magnetic  point  current  sources  in 
unbounded  anisotropic  uniaxial  media  is  presented.  The  Green’s  function  consists  of  an 
integral  over  the  plane-wave  spectrum  of  a  continuum  of  eigenvectors,  in  which  the  eigen¬ 
vectors  are  guided  (propagated)  in  the  preferred  z-coordinate  direction.  This  development 
uses  the  Lorentz  reciprocity  relation.  The  dyadic  Green’s  function  for  anisotropic  uniaxial 
multi-layered  media,  3m>0(r,  r'),  will  be  derived  in  section  2.4,  using  again  the  reciprocity 
relation  in  the  multi-layered  media.  In  Section  2.5  the  Green’s  dyadic  developed  in  the 
previous  section  will  be  used  to  derived  the  dyadic  Green’s  function  for  the  periodic  electric 
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as  well  as  magnetic  point  currents  in  anisotropic  uniaxial  multi-layered  media. 


2.2-A  General  dyadic  Green’s  function  formulation  for  anisotropic  media 

Maxwell’s  equations  in  a  source-excited,  unbounded,  anisotropic,  medium  with  (eJut)  un¬ 
derstood  and  suppressed,  can  be  written  in  the  operator  form 


£r(r,r')  =  -J 

where  the  dyadic  operators  are  defined  as 


(2.1) 


C 


(  jueQe  -V  x  I  \  /e\ 

V  Vxl  ju» oA  )  ’  V  H  )  ' 


(2.2) 


and  e,  p,  are  respectively  the  relative  permittivity  and  permeability  tensors  (with  respect 
to  free  space)  of  the  medium,  I  is  a  unit  dyad,  and  J  represents  the  impressed  electric 
as  well  as  magnetic  current  sources.  The  electromagnetic  field  column-vector  £  can  be 
represented,  using  the  electromagnetic  dyadic  Green’s  function  theory,  as 

r(r)=  J  dr'  E(t,t')  •  «7(r')  ,  (2.3) 

where  the  electromagnetic  dyadic  Green’s  function,  E( r,  r'),  is  2  X  2  vector  matrix  of  electric 
and  magnetic  dyadic  Green’s  functions,  and  v  contains  the  source  region.  In  this  general 
formulation  the  electric  as  well  as  magnetic  dyadic  Green’s  functions  due  to  the  impressed 
electric,  e,  and  magnetic,  m,  current  dipole  sources  are  both  considered.  If  J  is  an  arbi¬ 
trarily  oriented  point  column  current  dipole  source, 

J  --  Vb{t  -  r') ;  V  =  (  Pe’f  +  Pe’*  ]  ,  (2.4) 

y  Pm,i  i  Pm,:  J 

where  t,  z  stand  for  transverse  and  longitudinal  directions  respectively,  then  the  electro¬ 
magnetic  field  may  be  viewed  as  a  distribution;  namely, 


27(r)  =  S(r,r')-7>  ,  (2.5) 

where  in  vector-matrix  notation,  the  above  equation  can  be  identified  as, 

(  <?ee(r,r')  ae"*(r,r') 

\  ame(r,r')  0mm(r,  r') 

The  first  and  second  superscripts  denote  the  nature  of  the  Green’s  dyadics,  (i.e.,  electric  or 
magnetic),  and  the  type  of  the  current  sources,  respectively. 
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Having  established  the  relation  between  the  electromagnetic  field  and  its  associated 
Green’s  dyadic,  one  may  construct  such  a  dyadic  Green’s  function  to  find  the  electro¬ 
magnetic  field  due  to  arbitrary  electric  and  magnetic  current  distributions  in  unbounded 
anisotropic  medium.  One  can  utilize  the  guided  wave  mode  theory  [6]  to  construct  such 
a  dyadic  Green’s  function  from  the  spectrum  of  the  guided  (propagated)  wave  in  an  arbi¬ 
trary  z-direction.  In  the  following  the  property  of  the  guided  wave  modes  for  anisotropic 
media  will  be  explained,  and  then  the  theory  will  be  applied  to  find  the  Green’s  dyadic  for 
unbounded,  as  well  as  multi-layered,  uniaxial  anisotropic  media. 


2.2-B  The  z-propagating  mode  description  in  anisotropic  media 


In  order  to  effect  the  z-propagation  representation  of  the  eigenmodes,  one  can  decompose 
the  field  operator  into  transverse,  Ve,  and  longitudinal,  Jj,  components,  (V  =  V(  +  ^), 
which  yields, 


(2.7) 


where 

=  =  /  jVt  X  I  \  rir  =  (  0  -z  X  I  \ 

V  -J'V‘  x  1  ^oA  )  ’  \  z  x  I  0  ) 


(2.8) 


Also  the  dyadic  I,  and  /x,  may  be  represented  for  a  general  anisotropic  material  as 

e  =  €tt  +  zelt  +  etzz  +  zzet  ;  jj.  =  jitl  +  zfitt  +  jitzz  +  zifiz  ,  (2.9) 

where  for  (0  =  e  or  /z): 

fitt  =  ^  ^XV  j  !  +  y 0VZ  ;  $zt  =  x0zx  +  y 0zy  ,  (2.10) 

After  decomposing  the  field  components  into  transverse  and  longitudinal  (i.e.,  E  =  E*  + 
z  Ez  ;  H  =  Ht  4-  z  Hz),  and  taking  a  dot  product  of  (2.1)  with  z,  in  view  of  (2.4),  (2.7)  and 
(2.9),  one  will  have 


Ez  =  -fe„.Et  +  — - V(.(ixH()-^«(r-r')  , 

(z  ~JU€0  (z 

(2-11) 

Hz  =  —  +  — - - V<  •  (Et  x  z)  — _EuL£ — 6(r  —  r')  . 

Hz  -JUHoHz  J^'HoHz 

Also,  upon  eliminating  the  longitudinal  fields  from  (2.1),  one  can  find  the  equation  for  the 
transverse  field  vectors  [6] 
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u>e0 


iu-  -^l-i-zxvA^^x 


oj2(qHo 


Hz 


Et 


~3 


r— Z  X  Vt  -  Z  X  Vt^ 

fo  A1*  J 


h‘+4*xh‘ 


■  r — —  +  4^1>«£(r  -  r')  -  •-1-  -z  x  VtPm2S(r  -  r') 


An  •  + 


J«o 


1  -  _  1 .  _ 

-z  x  Vt— z  •  Vtx 


U/IqUz 


(2.12) 


(jJ^CqUcXo 


( Z 


+3 


^z  x  Vt  -  z  x  V£— 1  .  Et  -  j~ 

Uz  € 0 J  OZ 


z  x  Et 


-Pm^(r-  r,)  +  ^-Pmz6( r  -  r')  +  — ?— z  x  VlPez6( r  -  r')  . 
J  3  Hz  ufO  fz 


If  the  anisotropy  of  the  medium  is  transverse  to  z;  €ez,£*t,  Htzi  a-11^  Hzt  identically  vanish 
and  (2.11)  and  (2.12)  will  be  simplified  to: 


Ez  =  — ^ - V(  •  (z  x  Ht)  — ~Z—6(t  -  r')  , 

-jwe jwe0e* 

Hz  =  — - - V,  •  (Et  x  z)  -  J&£- 6(t  -  rO  , 

-3uPoHz  3ujHoHz 


(2.13) 


OJfo 


zxetf  +— 5 - V,  (-Z  •  Vtx  ) 


u}*eoPoHz 


E,-jrB,  = 


^P-^(r-^  +  -l-V,P„<(r-r,)  , 

3  u>PoHz 

(2.14) 


W/i0 


X  1  ^ 

z  X  Hu  ‘  +  .2.  T7  Vf(-z- Vtx  )|  Hf  -  j— Et  = 


u‘(oHo(z 


dz 

z  x  pmt«5(r  -  r')  1 


-V,pmz<5(r-  r')  . 


To  obtain  the  mode  function  characteristics  of  the  guided  waves  propagating  along  z 
for  the  (transversely)  unbounded  medium,  one  will  seek  the  source-free  solutions,  ('P  =  0), 
of  (2.1)  with  the  generic  form  of 


na(r,kt)  =  #a(p)e->*°‘ 


*«(P)  =  ^*ae~jk,  p  , 


(2.15) 
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where  'i'a(p)  and  Ka  are  the  solenoidal  eigenvectors  and  eigenvalues  of  the  operator  (2.7), 
respectively,  and  kt  is  the  transverse  vector  wavenumber  (i.e.  kt  =  >/kt  •  kt). 

For  the  medium  with  a  reflection-symmetric  property  (with  respect  to  a  plane  transverse 
to  z)  [6],  (i.e.,  the  anisotropy  confined  in  the  transverse  plane),  if  '&a>(p)  is  the  eigenvector 
corresponding  to  then  ^a<(p)  =  X?P,Q>(p)  is  the  eigenvector  for  -kq,  where  the 
reflection  operator  it  is  defined  as, 


X 


0  ) 
-lt  +  zz  J 


(2.16) 


where  I*  is  the  transverse  unit  dyad.  This  can  be  seen  from  (2.7)  and  (2.16)  (in  view  of 

=  =======  > 

(2.13)  and  (2.14)),  where  TtK  =  XX  and  XW  =  — W it.  Hence,  the  eigenvectors  ’&a<(p) 

corresponding  to  ±k„  can  be  expressed  as, 


* a<(p ) 


2  x 


P 


(  Et,a  i 


(2.17) 


where  e„,ha  and  eIia,hIj0(  are  normalized  mode- vectors  transverse  and  parallel  to  z ,  re¬ 
spectively,  for  the  eigenvalues  of  ±/cQ,  and  Za  is  the  mode  impedance  [6,  14].  Thus,  the 
eigenvectors  !PQ>(p)  and  !Pra<(p)  for  ( a  —  1,2)  can  be  written  in  terms  of  identical  mode- 
vectors  which  do  not  distinguish  between  the  wavenumbers  na  and  -Ka. 


The  analytical  expression  for  the  wavenumber  na,  and  the  eigenvector  ^(p)  for  (a  = 
1,2),  for  a  given  transverse  vector  wavenumber  kt  and  the  angular  frequency  of  uj,  are 
directly  related  to  the  relative  permittivity  and  permeability  tensors  (e,p).  Explicit  ex¬ 
pressions  for  Ka,  ea,  ha,e*i0,  h*a,  and  Za  for  a  general  uniaxial  anisotropic  medium,  where 
(i  =  (tI t  +  (Ziz  ;  p.  =  ntlt  +  pzzi)  are  given  by,  (see  also  Appendix  A  for  the  derivations) 


ei 

kxx  +  kyy 
kt 

hi 

kx  y  -  ky  X 

-  zxei=  fct  , 

(2.18) 

^2,1 

=  — —  ;  hz, !  =  o  , 

KX 

=  JkaPtft  -  ~k 1  ;  Zx  =  , 

V  (z  ve0  (t 

(2.19) 

kx  y  -  ky  X 

kt 
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kx  x  +  ky  y 

kt 


(2.20) 


h2  =  z  x  e2  =  - 


e*,2 


«2 


0  ;  = - - —  , 

Jklmtt  -  —k2t  ;  Y2  =  i- 
V  /i*  ^2 


*2 


(2.21) 


where  Ajq  =  u2CoHo,  and  A:t  =  ^Jk2  4-  k2.  Note  that  a  =  1  and  2  axe  respectively  correspond 
to  transverse  magnetic  TM  (or  .E-mode),  and  transverse  electric  TE  (or  E-mode)  respec¬ 
tively.  In  the  following  section,  the  dyadic  electromagnetic  Green’s  function  for  a  general 
uni-axial  anisotropic  unbounded  medium  will  be  formulated,  and  subsequently  the  dyadic 
Green’s  function  for  a  uniaxial  anisotropic  multi-layered  medium  will  be  derived. 


2.3  Dyadic  Green’s  function  for  unbounded,  uniaxial  anisotropic  media 

In  this  section,  the  theory  of  the  z- propagating  (guided)  eigenmode  for  a  general  uniaxial 
anisotropic  medium,  which  has  been  established  in  the  previous  section,  will  be  applied 
to  obtain  an  explicit  expression  for  the  electromagnetic  dyadic  Green  function,  S(r,r'), 
which  is  associated  with  the  electric  and  magnetic  point  sources,  r  —  r'),  in  unbounded 
uniaxial  anisotropic  media.  In  the  the  next  section  the  procedure  developed  here  will  be 
extended  to  obtain  the  dyadic  Green’s  function  for  the  uniaxial  anisotropic  multilayered 
medium. 

The  geometry  of  the  problem  dealing  with  full  space  uniaxial  anisotropic  medium  with 
the  relative  constitutive  parameters  of  (I  =  etIt  +  €zii  ;  /t  =  //tI t  +  //Zzz)  excited  by 
VS(r  -  r')  is  shown  in  Figure  2.2.  In  the  following  formulation,  the  notation  <  means  that 
the  entire  space  consists  of  two  region,  z  >  z'  and  z  <  z’\  and  as  shown  in  Figure  2.2,  z  =  z' 
is  the  plane  5  (normal  to  the  z-axis)  containing  source,  VS( r  —  r').  The  z-propagating 
representation  of  the  dyadic  electromagnetic  Green’s  function  for  the  full  space  anisotropic 
medium  in  terms  of  the  solenoidal  vector- wave  functions  can  be  expressed  as  [6,  14], 

-(r,r')=  /  dkt  S<(kt,r,r')  +  zz0z  ;  0*  =  (  ?  ^(r-r')  (2.22) 

7  \  U  -iw/Xo/ix  / 

> 

where  the  electromagnetic  dyadic  vector- wave  function  5  (k*,r,  r')  is  defined  as 

>  2  >  > 

S<(k„r,r/)=^  i7<(kt,r)^(kt,r')  ,  (2.23) 

a=l 
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> 

and  once  again,  the  source-free  solenoidal  vector- wave  function  J7<(k*,  r)  can  be  expressed 
as 


/7<(kt,r)  =  -Lva<e-ik'P*iKa* 
2ir 


(2.24) 


Outside  the  source  region,  z<z',  the  Green’s  dyadic  can  be  represented  in  terms  of  the 

> 

spectral  integral  of  the  source-free  solenoidal  vector-wave  functions,  eigenvectors,  J7<  as  of 
(2.23).  In  order  to  make  the  representation  complete  in  the  source  region  at  r  =  r',  the 
correction  term,  0*,  can  simply  be  added  to  the  spectral  integral  of  the  solenoidal  vector- 
wave  representation  of  the  electromagnetic  dyadic  Green’s  function,  which  is,  of  course, 
valid  outside  the  source  region,  [13].  The  variable  kt  (i.e.,(-kkx  +  yky);dkt  =  dkxdky) 
are  the  continuous  transverse  eigenvalues  which  span  the  entire  spectral  domain  (— oo  < 

kx  <  oo;  — oo  <  ky  <  oo),  and  the  summation  is  over  the  discrete  eigenmodes  (a  =  1  and 

> 

2)  propagating  in  the  z  direction.  The  unknown  vector  spectral  amplitude  «4^(kt  ,  r')  is 
defined  as  a  row  vector, 


(2.25) 


where  the  superscript  e  and  m  stand  for  electric  and  magnetic  sources,  respectively,  and  will 

be  found  from  the  procedure  using  the  reciprocity  relation  in  uniaxial  anisotropic  media, 

>  > 

(see  Appendix  B),  to  the  fields,  E  =  3  •  V,  of  (2.23)  and  the  source  free  solenoidal 

< 

vector-wave  function  of  I7>(— kj,  r)  in  the  region  V,  bounded  by  planar  surface  S>  and  5<, 
which  are  slightly  above  and  below  the  surface  5  of  the  Figure  2.2.1  respectively,  [22,  13, 
pages  2U0-206]1.  Hence, 


js  *^J7>(-k„r),  Wi7^  -jf  ds(n>(-kt,r),  W  27^  = 

-  Jv  dv  (vS(t-t'),  U  (n>(-kt,r)yj 


(2.26) 


=  -\V,U  JT>(-kt,r') 


where  U  (IIa)  is  a  diagonal  vector  matrix  comprised  of  the  elements  of  the  vector- wave 
mode,  IIa,  (see  (2.15)  and  (2.17);  tJ ( Tla )  and  W  axe  defined  as, 

;  W  = 

‘The  reciprocity  theorem  is  applied  to  the  volume  V  here  with  the  radiation  condition  implied  as  p  — >  oo. 


(2.27) 


^(iT«»=(o  -°h) 
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and  the  operator  notation  of  (II  a  ,  W  II $)  denotes  the  usual  scalar  product  of  the  two 
column  vectors: 


(J7«,  Wi7^)  =  iI^-'SVi7/3  =  Ea.(H/3xz)-H0.(zxE/3)  .  (2.28) 

Incorporating  (2.23)  into  (2.26)  and  using  the  orthogonality  relations  developed  in  Appendix 
B,  (see  (B.8)  and  (B.9))  yields, 


A<(kt,r') -V  =  (  a*  -pe 


1  VT.C(*a>(- kt))c+lkt.0,±iKn2, 
2t  -2  Za 


(2.29) 


Therefore,  from  (2.29),  (2.23),  (2.17)  and  (2.22),  the  2-propagating  plane-wave  spectrum 
(PWS)  integral  representation  for  the  full-space  anisotropic  uniaxial  dyadic  electromagnetic 
Green’s  function  can  be  identified  as. 


S(r,r')=J  ^5  T  +  «e,  . 


> 

r.  ^  • 


where,  the  vector-wave  of  is  identified  as, 


(2.30) 


*»(k') s  (» • ■ 6  + fe.  )r = ( f  •  «2-3i) 


and  1  is  a  dyadic  unit  column  vector- wave,  (i.e.,  1  •  E  =  E),  and  T  stands  for  transpose. 
The  PWS  representation  of  the  Green  ’s  dyadic  for  uniaxial  anisotropic  media  resembles 
in  many  ways  its  counterpart  of  the  dyadic  Green’s  function  for  isotropic  media  [13].  The 
Green’s  dyadic  S(r,  r')  can  be  decomposed  into  E-  and  H -  type  modes  with  the  different 
propagation  constants  (eigenvalues)  of  «Q(a  =  1,2);  moreover,  the  eigenmodes  satisfy  the 
reflection  symmetry  with  respect  to  a  plane  transverse  to  the  ^-direction. 


The  (PWS)  of  the  dyadic  Green’s  function  S( r,  r')  given  above  provides  information 
on  the  general  form  of  the  (PWS)  solution  for  the  electromagnetic  fields  for  multi-layered 
uniaxial  anisotropic  media.  In  the  following  the  dyadic  Green’s  function  for  the  uniaxial 
anisotropic  slab  will  be  derived,  which  then  will  be  generalized  to  multi-layered  uniaxial 
anisotropic  media. 


2.4  Electromagnetic  dyadic  Green’s  functions  for  uniaxial  anisotropic 
multi-layered  media 

The  presence  of  a  planar  boundary  in  anisotropic  media  generally  introduces  coupling  be¬ 
tween  the  wave-vector  modes  which  are  traveling  from  and  toward  the  source,  and  hence 
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makes  any  one  amplitude  dependent  on  the  amplitudes  of  the  others,  [14].  Nevertheless, 
for  uniaxial  anisotropic  multi-layered  media  the  mode  couplings  do  not  occur. 

The  electromagnetic  dyadic  Green’s  function  for  uni-axial  anisotropic  multi-layered  me¬ 
dia  can  be  expressed  as  sum  of  E(r,  r')  in  (2.30)  and  another  contribution  to  account  for  the 
field  scattered  by  the  layered  media.  The  scattered  field  contribution  can  be  expressed  in 
terms  of  a  PWS  integral  resembling  that  for  S( r,  r').  Consider  an  arbitrarily  oriented  point 
column  current  dipole  source  as  of  (2.4)  located  in  region  (n)  with  constitutive  parameters 
(ln  =  fi.nlt  +  cz  nzz  ;  p.n  =  Ht,nh  +  Hz,niz)-  Let  us  represent  the  Green’s  dyadic  for  region 
m  as 


/_  > 

dkt  —  <m(kt,r,  r')  +  zzO*,n  , 

where  the  spectral  Green’s  dyadic  within  the  region  m  can  be  identified  as 


(2.32) 


> 


*<,n(kt,r,r')  =  £17? >(lct,r)X)<r<(kl,r')  , 


(2.33) 


a=l 


with  n„  is  defined  as 


J7?V<,r)  =  #m<(kf,z)e--’kt*>  ;  #m<(k,,  z)  =  ^*«m<e**‘?*  +  R?  < (0)¥am> e±jK? ^  , 

(2.34) 

> 

where  subscripts  (m,n)  denote  the  field  and  source  regions,  respectively;  R™<(z)  is  the 
effective  modal  reflection  coefficient  at  the  interfaces  (m,m  +  1)  and  (m,m  -  1)  for  (>) 
and  (<),  respectively  evaluated  at  z,  and  the  superscripts  <  denote  that  the  incident  field 

travels  in  ±  2-directions;  hence,  (i.e.,  R™  <{z,)  =  R™  <(2)t)  e±j2*™ (*)-*'•'>).  The  unknown 

> 

V is  the  vector-wave  dyadic  amplitude,  for  region  m,  and  is  defined  as  a  row  vector, 

> 


T>am<  =  (  d%  d-  ) 


> 

m< 


(2.35) 


where  subscripts  e,  and  m  denote  the  electric  and  magnetic  type  sources  respectively,  and 
will  be  determined  from  the  Lorentz  reciprocity  relation  shortly.  The  subscript  <  here 

>  > 

denotes  the  region  of  z<z' .  The  effective  reflection  coefficient  R™<(z)  for  region  m,  is  a 
function  of  all  successive  layers,  (i.e.,  m  ±  l,m  ±  2, ^  ^  ^ )  of  the  uni-axial 

anisotropic  multi-layered  media,  and  is  given  by  the  following  recursive  relation  [13] 
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(2.36) 


>,  _  T?<  +  R?±'<(zm±l)e-^±l<1^ 

\zrn)  —  > 


1  +  IT  <  R?±x  <(zm± 


where  the  reflection  coefficient  r™<  is  given  by 


<  =  Z  ~  Za 

Za±l  +  Z™ 


(2.37) 


and  the  mode-impedances  Za  are  defined  (see  also  (2.18),  (2.20),  and  Appendix  A)  as 


rrn  *1  .  >7n 

6\  —  -  ,  £ 2  ~  n  - 

W  (0(t,n  «2 


(2.38) 


>  > 

Also  the  modal  dyadic  coefficients  13  ”  <  of  region  m,  and  <  of  region  n  on  either  side 

> 

of  the  source  are  related  via  the  effective  transmission  coefficient,  T<  n.Q,  and  is  given  by 

[13]2 

>  mil  > 

T<n;a  =  n  ra‘  <(-Ji)  e~}K  *  ;for  n<m  ,  (2.39) 

i=n 

> 

where  the  transmission  coefficient  Tm<  at  the  interface  (m,  m  ±  1)  is  given  by 


r„”  <(*„)  = 


i  +  r”< 

> 


1  +  T™  <  R?±l  <  *»±i 


(2.40) 


Hence,  one  only  needs  to  find  the  modal  coefficients  'D<  in  region  n,  in  order  to  com¬ 
pletely  specify  the  fields  in  all  regions.  Specifying  (2.33)  for  region  n  (i.e.,  m  —  n)  and 

>  > 

invoking  the  Lorentz  reciprocity  theorem  to  the  fields,  Sn  n  —  Sn  n  •  V,  of  (2.33)  and  the 

< 

source  free  solenoidal  vector-wave  function  of  II a  (-kf,  r),  of  (2.34),  in  region  n  of  volume 
V,  bounded  by  planar  surface  S>  and  5<,  which  are  slightly  above  and  below  the  surface 
S  of  the  Figure  2.1.2  respectively,  yields 


JNote  that  T£no  is  the  ratio  of  the  modal  incident  wave-amplitudes  of  regions  m  and  n;  namely, 


p  —  t<  f>n  <  P 

r  —  L  m.n.a  r  - 
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w ds  (nna>{-*t,r),^  i£n) 

-JvdV{ 


V6(r-r'),U  \n”(-kt,r) 


(2.41) 


Incorporating  (2.33)  into  (2.41)  and  using  the  orthogonality  relations  developed  in  Appendix 
B,  (see  (B.8)  and  (B.9))  yields, 


^"V„r') 


=  U 


Pe  Ac 


r = 


i  vT 

2t r 


U(*n>(-kt,z')):+l]ti.0, 


-2  AS 


(2.42) 


where 

AS  =  ZS  (1  -  R^iz^R^izn^e-2^^)  .  (2.43) 

Therefore,  from  (2.42),  (2.33),  (2.17),  and  (2.31)  the  z- propagating  plane-wave  spectrum 
(PWS)  integral  representation  for  the  electromagnetic  dyadic  Green’s  function  of  (2.32)  can 
be  identified  as. 


=•  (r  r')-  /  ^  t  T<  *«mW)»r(k e 

—  m.ntr,  T)~  J  4jr2  Am,n;o  -2  A"  +  ZZU*n  • 


(2.44) 


where  once  more  the  wave-vectors  #  and  #  are  given  by 


Va<e*jK°x  +  R<( 0)*a>e±jK°*  , 


> 


>  >  < 

+  RZ(0)#ae±jKaZ'  . 


(2.45) 


The  resemblance  of  the  electromagnetic  Green’s  dyadic  for  the  uni-axial  anisotropic  multi¬ 
layered  media  of  (2.44)  to  its  associated  counterpart  for  the  full-space  uni-axial  anisotropic 
medium  should  be  noted.  Both  dyadic  Green’s  functions  are  comprised  of  the  vector-wave 
mode  functions,  which  are  the  solutions  of  the  source-free  fields  in  the  uni-axial  anisotropic 
medium.  In  the  full-space  case,  the  excited  modes  are  those  which  travel  away  from  the 
source,  whereas  in  the  multi-layered  case  the  modes  which  travel  toward,  as  well  as  away 
from,  the  source  will  be  excited. 
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In  the  following  explicit  expressions  for  the  electric-electric  Qee ,  magnetic-electric  Qme , 
electric-magnetic  Gem ,  and  magnetic-magnetic  Gmm  dyadic  Green's  functions  will  be  presented3. 
In  view  of  (2.17)  and  (2.31)  vector-wave  functions  i>  and  #  of  (2.45)  can  be  written  in  terms 
of  the  electric  and  magnetic  mode  functions  as 


> 

#<  = 


(Zaea  ±  e*iQ )  e*iKat  +  R<(0)  ( Zaea  ^  ez,a)e±>K°z 
(±hQ  +  +  R<( 0)  (ThQ  +  hz,a)  e±JKaZ 


(2.46) 


(Zaea  t  e,^)e^«“*  +  R<(0)  ( ZQea  ±  ez,a) 

(Th a  +  hz,Q)e^  +  7Z<(0)(±ha+  hz,a)e±^i 


(2.47) 


therefore  the  electromagnetic  dyadic  Green’s  function  can  be  identified  as 


,p  pg  _  (  rO  Gem(ry)  \ 

— m.nir,  r )  -  ^  0«.(r>pf)  Gmm(r,r')Ji 


(2.48) 


I  4^2  « _ o  An  ^  ^ 

•  v  (kt.**)  #^<(k„*)*:,;(k*.*#) 


+  zz  (  -J40£i  n  ?  )  «( 

V  -J*HO  Ml.n  / 


r-r')  . 


CTe  /  A  _  f  -jk ,  (p-p')  Y' 

Vmnl  ’  ^  “  J  4ff2  ^  -2 A! 


,  -2  AS 

a=l  a 


(2.49) 


((Z°ea  ±  6,’° 


)m  2  +  JC<(0)  (Z0e«  =F  eXt0)m  e**™  2  • 


(Zaea±ez,a)ne±’*°*'  +  J£>(0)(ZQea 


+  zz - ^(r-r')  . 

-jwc  0elt„ 

’Note  that  the  expression  in  (2.44)  is  a  compact  representation  of  the  electric  as  well  as  magnetic  dyadic 
Green's  function  due  to  the  electric  and  magnetic  point  sources  in  a  uni-axial  anisotropic  multi-layered 
media. 
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<C!„(r,r')  = 


fQSi  -jk ,  (P-P')  y' 

4tt2  "  -2  A" 

Qr=l  a 


(2.50) 


(iha  +  +  /C<(  0)(Tho  +  hz,a)me±2"«* 


(Zaea  ±  et(Q)n  «***''  +  J£>(0)  (ZQeQ  T  ez,a)n  e*^*' 


GVr.n(r,r')  = 


fOSi  -jk,  (P-P')  V- 

4?r2  ^  -2  A2 

Q  =  1  a 


(2.51) 


(ZaeQ  ±  ez,a )m e*><*  +  ^<(0)  (ZQea  T  ez,Q)m  • 


(±ha  +  h^e^*'  +  iZS>(0)(ThQ  +  hz,a)ne*2< 


^^(r,r')  = 


f  fOfi  e-jMP-P')  V'  ZlLuif 
J  4x2  “  -2A{ 

Q  =  1  * 

|(±ha  +  h,a)ra€^l+. 


,  ~2  AS 
0  =  1  a 


(2.52) 


,)m  +  C<(0)  (Tha  +  hz,a)m  e*2"™  *  • 


(±hQ  +  hz,Q)„  e^<z  +  ^>(0)  (Tha  +  hz,a)n  6**2* 


4.  iz — : - £(r  -  r')  . 

-JvHoHz,n 


2.5  Electromagnetic  dyadic  Green’s  functions  for  uniaxial  anisotropic 
multi-layered  media  due  to  periodic  electric  and  magnetic  point  cur¬ 
rents 


In  the  following  the  dyadic  Green’s  function  for  an  infinite  periodic  array  of  an  arbitrarily 
oriented  current  point  source  in  a  uni-axial  anisotropic  multi-layered  media  will  be  derived. 
The  geometry  of  the  problem  is  shown  in  Figure  2.3. 

The  starting  point  for  this  analysis  is  based  on  the  following  theorem.  Consider  a 
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periodic  impulse  train,  /  with  an  spatial  phase  variation  uq  and  period  X : 


/(*)  =  Yi  6(x  -  nX)e~J“0X  ■  (2.53) 

n=  —  oo 

the  Fourier  transform  of  this  function, denoted  by  F  is  given  by 

F(u>)  =  J  dx 

OO 

=  E 

n=  — oo 
oo 

-  E 

n=— oo 

but  from  the  distribution  theory  [16] 

E  e-"T  =  Y  t  ,  (2-56) 

n=— oo  n=— oo 

so,  it  follows  that, 

n«)  =  |  f;  «((w-«o)-?p) , 

n=  — oo 

which  simply  state  that  a  periodic  impulse  train  has  a  Fourier  transform 
periodic  impulse  train. 

Let  us  assume  that  the  current  elements  are  spatially  periodic  with  the  period  of  X' 
and  Y'  in  x'  and  y'  directions  respectively,  where  {x',y')  coordinate  are  skwed  a,  and  (3 
degrees  with  respect  to  ( x ,  y)  global  coordinate  respectively,  as  shown  in  Figure  2.3.  Thus 
the  location  of  (i,g)th  element  OP  can  be  represented  as 


(2.57) 
which  is  also  a 


e’“xf(x) 


J  dx  ej{“-“°)x6(x-nX) 


j(u-u0)nX 


(2.54) 

(2.55) 


OP  =  iX'x'  +  qY’y' 

=  i  X'(xcosa  -  y  sin  a)  +  ? ^(xsin /?  +  y  cos/9)  (2.58) 

=  (i X'cosa  4-  qY'sin(3)x  +  (— * X'sina  -f  qY' cos(3)y  . 

Let  also  assume  that  the  current  elements  are  phased  at  an  angle  ( 9,4> );  thus,  the  discrete 
periodic  current  distribution  can  be  denoted  by 

OO 

J(  r,r')  =  E)  Ve-}U*X-J“'>y6(z  -  z')  • 

«,<?=-oo 

f)  (x  -  (i  X1  cos  a  +  qY1  sin  (3))  6  (y  -  (-iX'sina  +  qY1  cos/3))  ,  (2.59) 
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where 


(2.67) 


In  the  light  of  the  above  discussion,  the  electromagnetic  dyadic  Green’s  function  for  the 
periodic  current  point  source  can  be  expressed  as, 


-p  _  (  Oee(r,r’)  0en*(r,  r') 

-m.nir.r;  ^  £me(r,r')  0mm(r,  r') 


=  J—  y 

X'Y'  .  ^ 


D-jk;«(P-P')  m 


(2.68) 


V  r^n;a(Kq)  (  #£<(1$*,*)' Ki  iK* ,  *')  *^0$*  i  *)  KM*  ’  *')  |  , 

2-i  1  >  <  >  <  .  i 


-jurf0f*.n 

0 


.OAn/lf'Wl  I  >  „  ^  I 

Q=1  a(  ‘  }  \  *ommW,*)<e  (k;* /)  *z<(K(l>z)*ZX(Kq^')  / 

o  \  °° 

!  j  6(z  —  z')  ■  y^  6  (x  -  (i'X'cosq  +  qY'  sin/3))  8  (y  -  (-iX'sina  +  qY'cosp))  , 

—jv[l  QHz,n  /  i,q=  —  oo 


where  subscript  p  denotes  the  periodic  type  Green’s  function,  and  k},,J  is  discrete  transverse 
eigen- value,  and  is  given  by  defines  as 


k |*  =  *i*x  +  *j*y  , 


(2.69) 


where  kl*  and  k'^q  are  defined  in  (2.66).  Likewise  explicit  expressions  for  the  periodic 
electric-electric  Gp:ee ,  magnetic-electric  £  p:  me ,  electric-magnetic  Cp:em,  and  magnetic- 
magnetic  G  p  mm  dyadic  Green’s  functions  can  be  expressed  as 


e:,e„e(r,r')  =  —  y 


,,-jkj  "  (P-P’)  y  ^m,n;o(^t<?) 

0=1  -2K(Kq) 


(Zaea  ±  eX'a)m  +  fC<(0)  (Zaea  *  ez,a)m  e**?* 


(ZaeQ  ±  eZ  Q)n  e±JK^'  +  ^>(0)  (Zaea  T  e,,a)„  e=^*'  + 


1 

zz — : - 8(z  -  z')  •  V'  £  (z  -  (i  A^cosa  +  gy'sin/?))  6  (y  -  (-i  X' sin  a  +  qY' cos  $)) 

-JWfO^z.n 


<?^:.ne(r,r')  =  -i_  ^  e-jk;’(p-p')  £ 


T<(„;o(k;’9) 


,^oc  0=1  -2AS(kJ’<I) 

(±ha  4-  +  iZ^<(0)(Tho  +  hz,0)rae« 


(2.71) 
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(Zaea  ±  e,,a)n  e±3K°z'  +  J^>(0)  (Zaea  T  e,,a)„  e*^*' 


Cpem(r  r')  =  — - —  V  e 

vm,n  V1)1  )  \'Y'  ^ 


c-jk \  q  (p-p‘)  ST"  '^‘m,n;a(^ct <?) 

«i  -2  AS(k^)  ‘ 


(2.72) 


i,q=-oo 

1  >  \ ,’9 
(Zaea  ±  e,,a)m  e^“-  *  +  i2™<(0)  (Zaea  T  ez,0)m  • 

f(±hQ  +  4-  R2>( 0)(^hQ  +  hfl0)neW"S*')  9 


(JP:mmfr  r')  -  _ _ 

y  m,n  V 1  i  1  )  Yt 


1 


oo 

X'K'  .  ^ 

1,0=  — OO 


-2As(k;'9) 


f(±hQ  +  hIiC 


)m  +  iC^(O)  (Th„  +  Va)m 


(■ 


'.<7 


(±ha  4-  ht,0)n  e^S*  +  J£>(0)  (ThQ  4-  h*,a)n  e:fJ*°2'  + 


1 


(2.73) 


zz__ - ^  6  (x  -  (iX'cosa  +  qY'smfi))  6  (y  -  (-iX'sina  +  qY' cos /3)) 


-JW/ioMz,n 


i,q—  —  oo 
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=e>  ~  (i,th  +  €iiX zz  ;  -f  (ii}Zzz 


Figure  2.1:  Electric  and  magnetic  point  current  dipole  sources  in 
anisotropic  uniaxial  multi-layered  media,  ( =e  =  etIt  +  ezii  ;  = 

/ztIt  +  fJ-ziz).  Also  the  planar  surface  S>  and  5<  slightly  above  and 
below  the  §ource  are  shown. 


Figure  2.2:  Imaginary  plane  S,  parallel  to  xy  plane,  passing  through  the 
source  at  z  =  z!  in  unbounded  anisotropic  uniaxial  media.  Also  plotted 
are  planar  surfaces  S>  and  5<  slightly  above  and  below  the  source. 


\  +  €i,zii  ;  A4.'  =  ZZ 


Figure  2.3:  Infinite  periodic  array  of  arbitrarily  oriented  current  point 
sources  in  a  uni- axial  anisotropic  multi-layered  media.  The  current 
elements  are  spatially  periodic  with  the  period  of  X'  and  Y'  in  x'  and 
y'  directions  respectively,  where  y')  coordinate  are  skwed  a,  and  {3 
degrees  with  respect  to  (x,y)  global  coordinate  respectively.  OP  is  the 
spatial  position  of  (i,g)th  element.  Also  current  elements  are  phased 
at  an  angle  ( 9,<f> )  (not  shown). 
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Chapter  3 


Electric  dyadic  green’s  function  for  a  grounded  uniaxial  anisotropic 
double  layered  slab  due  to  an  electric  point  source 


3.1  Formulation  of  the  Sommerfeld  representation  for  the  dyadic  Green’s 
function 

The  planewave  spectrum  (PWS)  integral  representation  of  the  electromagnetic  dyadic  Green’s 
function  for  a  general  uniaxial  anisotropic  multi-layered  medium,  which  has  been  derived 
in  the  previous  chapter  will  be  specialized  to  a  grounded  double  layered  unaxial  anisotropic 
slab  as  shown  in  Figure  3.1,  where  an  electric  point  source,  (pe<5(r  -  r'),  pe  =  pfe  +  pie), 
is  located  in  region  (2).  Then  the  Sommerfeld  2-propagation  integral  representation,  and 
subsequently,  the  radially  propagating  integral  representations  of  the  electric  dyadic  Green’s 
function  will  be  derived.  The  Sommerfeld  integral  representation  leads  to  the  singularity 
removal  technique  which  is  very  efficient  when  the  observation  point  is  very  close  to  source 
point.  Nevertheless,  as  the  lateral  separation  between  the  source  and  observation  points 
increases,  the  Sommerfeld  integrand  starts  to  oscillate  rapidly,  and  eventually  the  method 
loses  its  efficiency.  These  oscillations  are  due  to  the  rapid  growth  of  the  argument  of  the 
Bessel  function  in  the  integrand  which  causes  the  singularity  extraction  method  to  eventu¬ 
ally  lose  its  efficiency.  The  radial  propagation  representation  of  the  dyadic  Green’s  function 
is  very  efficient  when  the  source  and  observation  points  are  laterally  separated.  This  ef¬ 
ficiency  is  due  to  the  imaginary  argument  of  the  Hankel  function  which  is  present  in  the 
radial  propagation  of  the  dyadic  Green’s  function. 

The  electric  field  is  related  to  the  electric  dyadic  Green’s  function  as 

Em(r)  =  J  Je(r')  ;  m  =  1  or  2  ,  (3.1) 

where  Je  is  the  impressed  electric  current  within  the  slab,  and  V  contains  the  source  region. 
If  the  J(r')  is  an  arbitrarily-oriented  point  current  source  of  pe6(r  -  r'),  the  electric  field 
may  be  viewed  as  the  distribution 

Em(r)  =  0"2(r,r')  •  pe  ,  (3.2) 

The  electric  dyadic  Green’s  function  due  to  the  electric  source  can  be  expressed,  using 
(2.49),  as 
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^m,2(r>r')=/  T~~J  e  jkt'(P  P)^m,2( kt,  Z,  Z')  +  ZZ——^ - 6(r  -  T1) 

J  ***  —]U€q€Z'2 


where 


^  *  ( P  P  )  —  kx(x  x  )  -f*  ky(y  y  )  J  d  k$  —  d  kx  d  ky  , 


and  Gm  2(k*,  z,  z ')  can  be  represented  in  the  matrix  form 


C<  (k  ~Gn(kt,z,x')  ~Gfx(kt,z,z') 

Vm,2(kt,Z,Z  )  -  I  >  > 


\  G;t(kt,z,z’)  G;z(kt,z,z’) ;  m  2 

the  first  and  second  subscripts  denote  the  electric  field  and  electric  point  current  directions, 
respectively.  Using  analytical  expressions  for  the  modal  fields,  ea,elQ,  ha,  hz0[  and  the 
modal  impedance  Za  given  in  Appendix  A,  the  Green’s  dyadics  inside  the  bracket  for 
m  =  2  (i.e.,  field  points  located  in  region  2)  are  explicitly  defined  as, 


Gft(kt,z.z') 


-4—  ^4;<(kt, z, z')k;k;  +  ki^g';<iku  ,  (3.6) 

—zu/e 0  ( t  «o 


where 


Gfz{kuz,z')  =  g<(kt,z,z')k'tz  , 

G*t(kt,z,z')  =  3u(kt,z,z’)zk't  , 

>\  kt  1  ,w 

=  ~2^7z-729  iKz'z)zz  ’ 


O  _  Xfcx  +  y  ky  f  „  _  ±  f,  y  kx  -Xky 

—  i  ,  Kt  —  Z  X  Kt  —  “ 


(3.10) 


and  the  functions  <7'(kt,  z,  z')  and  g"(kt ,  z,  z')  are  defined  in  Appendix  C.  The  prime  (/)  and 
double  prime  (//)  here  are  respectively  correspond  to  a  =  1,  and  a  =  2. 


The  Green’s  dyadic  cam  be  written  in  terms  of  scalar  functions  if  one  notes 


kj  ’  ‘  * 


(z  x  Vt)(z  x  Vt) 

kt 


(3.11) 


and  uses  the  identity  of 
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v,v,  —  <*£+*£x*3j+y£) 

..a2  ..a2  ...  ...  a2 

=  xx^  +  yyj^  +  (xy  +  yx)^  , 


(3.12) 


(z  X  V,)(z  X  Vi)  —  xiJL  +  yyJL-fry+y*)^- 

a2 .  a2 .  ...  ...  a2 

s  xx<v'-^>  +  yy<v'- V,‘(xy  +  yx)a7a7 


(3.13) 


Thus,  the  planewave  spectral  dyadic  Green’s  function  (3.3)  (in  view  (3.5)-(3.9))  can  be 
expressed  as 

S,1(r,r')  =  ji-|xx^-tJ«<(r,r')  +  ^W<(r,r')j 


y(-k10U<(r,r')  +  -~Wi(T,r') 
+(xy  +  yx)^j*V<(r,r')}  , 


c“(r'r,)  =  5^7 Vl  pr’r')] 4  • 

[&'•*  • 
6>Uty)  =  -2^Skrv?(r’r,)  • 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


where 


g<(r,T>)  =  J  jk*^  pn>  , 

U<(r,r')  =  j  ^e-*‘  (p-p')^g"<(kt,z,z')  , 


(3.18) 


(3.19) 


30 


(3.20) 


W<(r,r/) 


V<(r,r') 


/  f!hi  e-J*.  (P  -  p') 

J  4*2 


-^~2 9u(kt,  z,  A  -  kl  z ,  z) 

(2  ,tk(  K2*Ct 


J 


c-jk,  (p  -  p') 

4*2  * 


The  Equation  (3.18)  has  a  Fourier- Bessel  representation  given  by 


(3.21) 


^(r’r')  =  H  £  e-J n(4>-4>')  £°  g<(t,z,z')Jn(kotp)Jm(koSp')Sd,t ,  (3.22) 


which  can  be  obtained  from  (3.18)  by  employing  the  following  transformations 


kt  =  k0£  ;  kx  =  ko£  cos  a  ;  ky  =  ko£  sin  a  (3.23) 

(x  -  x')  =  p  cos  <(>  -  p'  cos  4>  ;  y  -  y'  =  p  sin  <p  -  p'  sin  <p'  .  (3.24) 

If  the  coordinate  system  is  chosen  so  that  p'  =  0,  then  the  expression  (3.22)  reduces  (in 
view  of  Jo(0)  =  1,  and  Jm{ 0)  =  0  for  m  ^  0)  to 

fl'(r’r)  =  ^*  ij  .  (3.25) 

The  Fourier- Bessel  representation  of  the  U,  W,  and  V  are  obtained  using  the  same  procedure 
as  of  g  and  are  given  by: 


W<(r,rO 

W<(r,r') 

V<(r,rO 


IT  Jq  K 2 

|r  I" 

IT  J o  I  C2,t  k2 


MZkop)^- 


^  f30  K&((,z,z^J0((*0P)(d(  . 

Z7r  Jo  k2 


(3.26) 


(3.27) 

(3.28) 


Expressions  in  (3.25)-(3.28)  are  the  Sommerfeld  type  integral  representations  of  g,  U,  W, 
and  V  given  in  terms  of  planewave  spectral  (PWS)  integral  representation  given  in  (3.20)- 
(3.25).  It  should  be  noted  that  the  wavenumbers  k.\  and  k "  in  (3.25)-(3.28)  are  normalized 
with  respect  to  ko,  namely  by 


=  \Jni  ~  K  f2  5  =  y/nt  -  n"  £2  , 


(3.29) 
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where  n,  =  eiitpiit,  n-  =  ^,and 

Each  term  in  the  integrand  of  (3.27)  exhibits  a  spurious  pole  at  £  =  0.  This  leads 
to  an  integral  with  a  logarithmic  singularity  and  it  is  a  consequence  of  interchanging  the 
derivative  and  integration  operations.  Strictly  speaking,  this  interchange  is  not  allowed 
(because  of  the  convergence  of  difficulty  as  £  -+  0).  However  the  expressions  in  (3.27)  as  a 
whole  converge,  which  means  that  the  integrand  of 

((,!,*')]  M(kop)  7  (3.30) 

*2,1  k2  < 

is  regular  at  £  =  0.  This  can  be  shown  by  explicit  substitution  of  g'tt  and  g"t  and  taking 
the  limit  of  the  integrand  of  the  above  expression  as  £  — >  0.  So  in  a  distributional  sense  an 
interchange  of  the  operators  is  justified. 

The  convergence  behavior  of  these  functions  in  (3.25)-(3.28)  depends  strongly  upon  the 

value  of  | z  -  z'\  (vertical  separation  between  the  source  and  observation  points).  For  the 

cases  of  planar  structures,  such  as  microstrip  arrays,  the  factor  which  exhibits  an  exponential 

decay,  I,  goes  to  unity.  Thus,  as  £  — ►  oo,  the  Sommerfeld  type  integrands  decay 

algebraically  as  ~  /  rather  than  exponentially;  the  algebraic  decay  results  from  the 

\JV«SP 

asymptotic  behavior  of  J0(£  k0  p)  for  large  £  k0  p. 

It  is  a  well  known  fact  that  the  numerical  integration  of  oscillatory  functions  are  very 
time  consuming,  particularly  when  they  have  an  end  point  at  infinity.  To  overcome  this 
difficulty,  which  is  inherent  in  any  open  structure  involving  multilayered  media  one  may 
use  the  singularity  removal  technique.  This  technique  is  very  efficient  when  the  lateral 
separation  (k0p)  is  small  (i.e.,  if  kop  <  1).  It  deals  with  asymptotic  behavior  of  the  inte¬ 
grands,  which  in  turn  is  dictated  by  the  behavior  of  J0(£fcop)-  These  asymptotic  values  can 
be  integrated  in  closed  form;  therefore,  if  one  subtracts  these  asymptotic  values  from  the 
integrands,  then  the  resulting  integrands  are  relatively  smooth  and  fast  decaying  so  that 
they  can  be  integrated  efficiently. 

The  singularity  removal  technique  [19]  which  improves  the  convergence  of  the  conven¬ 
tional  Sommerfeld  integral  representation  is  described  as  follows. 

Let  I  be  a  typical  value  of  g,  U,  W,  or  V  of  (3.25)-(3.28)  which  in  general  has  the  form 

/=fr  F(Z,  z,  z')  J0{£  kop)  d  £  ,  (3.31) 

Jo 

where  £m  F(£,z,z')  is  recognizable  from  (3.25)  through  (3.28).  As  £  — *  oo,  the  value 
of  F(£,z,  z’)  asymptotically  approaches  some  complex  value,  say  5Z„Cne-^n;  (where  the 
explicit  expressions  associated  with  g,  U,  W,  and  V  can  be  found  from  (3.25)-(3.28),  and 

,>  „> 

the  limiting  values  of  g,f{£,  z,z')  and  gif{£,z,z')  as  £  — ►  oo  which  are  given  in  Appendix 
D.  But 
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(3.32) 


j»oo  Jrn  /*oo 

l  M(kop)e-^‘d( 

=  (-l)m — — — 

v  '  d(M)m 

Therefore 

/•OO  ,7™  /  1 

l~L  (W)-EC.. -*«)«■  *(«*.,) 

(3.33) 

Mathematically,  Equations  (3.31)  and  (3.33)  are  equivalent;  numerically  however,  the  in¬ 
tegral  of  (3.33)  possessed  a  much  smoother  function;  therefore  an  eight  or  sixteen  point 
Gaussian  integration  routine  is  sufficient  to  achieve  a  reasonable  accuracy  in  evaluating 
such  an  integral.  More  discussion  about  the  behavior  of  integrands  in  (3.31)  and  (3.33)  will 
be  presented  in  Chapter  5  which  is  devoted  to  numerical  results. 

It  should  be  pointed  out  however  that  the  smoothness  of  the  integrand  in  (3.33)  dete¬ 
riorates  as  the  lateral  separation  of  the  source  and  observation  popints  increases  (i.e.,  as 
k0p  becomes  large).  The  second  method  involves  the  use  of  a  radially  propagating  inte¬ 
gral  representation  wh;ch  overcomes  the  difficulty  associated  with  the  first  method,  and 
converges  very  fast  if  the  source  and  observation  points  are  separated  laterally  rather  than 
vertically.  This  technique  also  leads  to  the  uniform  asymptotic  closed  form  solution  which 
will  be  discussed  in  the  following  chapter. 

3.2  Alternative  radial  p-propagation  representation  for  the  dyadic  Green’s 
function 

The  starting  point  in  the  derivation  of  the  radial  propagation  microstrip  Green’s  function 
is  Equation  (3.25)  which  constitute  the  Sommerfeld  representation  of  g(r,  r')  as 

<7(r>r0  ~  2^  J0  Jo  Wop)  t  d  £  .  (3.34) 

The  above  integral  in  (3.34)  can  be  written  in  terms  of  the  Hankel  function  where  the 
integral  has  the  lower  limit  as  oce~]* ,  and  the  upper  limit  at  oo,  provided  that  g  is  an  even 
function  of  £.  Thus,  (3.34)  becomes 

5(r,r')  =  ^f  g(t,zt2')H<2\tk0p)tdt  ,  (3.35) 

JO O  C~~jW 


33 


in  which  the  contour  of  integration  is  such  that  the  branch  point  at  £  =  0  introduced  by 
the  transformation  is  avoided  as  shown  in  Figure  3.2. 

An  alternative  ^-propagation  representation  of  g  can  be  derived  by  deforming  the  inte¬ 
gration  contour  of  (3.35)  so  that  it  encloses  all  the  singularities  of  the  z-dependent  function 
g(£,  z' )  in  the  integrand.  General  techniques  for  arriving  at  alternative  representations  for 
Green’s  functions  have  been  discussed  in  [6,  18]. 

In  order  to  enclose  al'  .he  singularities  of  g(£,  z,  z')  in  the  f  plane,  the  analytic  property 

>  > 

of  g(£,  z,  z')  needs  to  be  studied.  As  mentioned  earlier,  explicit  expressions  for  g M<  and  g,<: 

> 

are  given  in  Appendix  C  (in  particular  see  (C.l)  -  (C.ll)  and  (C.32)  -  (C.39)).  The  g "<  and 

> 

g,<:  are  functions  of  the  propagation  constants  Kmi  m  =  0, 1  and  2,  and  in  the  complex 
£  plane  they  have  a  pair  of  branch  points  at  f  =  ±1  and  a  pair  of  branch  cuts  from  £  =  ±1 
to  infinity.  These  branch  cuts  must  be  defined  in  such  a  way  that  Im  (i/l  -  <  0  on  the 

entire  top  sheet.  This  is  the  only  pair  of  branch  cuts  that  allow  the  analytic  continuation 

>  > 

of  g"<  and  g ,<:  into  the  entire  top  sheet  of  £  such  that  they  satisfy  the  radiation  condition. 
Therefore,  the  deformation  of  the  real  axis  path  of  integration  into  the  complex  plane  is 
possible  as  long  as  it  stays  on  the  top  Riemann  sheet  as  |f|  — »  oo;  of  course,  one  must 

>  > 

account  for  residues  associated  with  any  poles  in  the  integrals  (g"<  and  gIK)  which  are 

enclosed  in  this  contour  deformation.  As  mentioned  previously,  Figure  3.2  shows  a  typical 

> 

complex  £  plane  containing  branch  cuts  due  to  the  Hankel  function  and  g"<  in  the  integrand 
of  (3.35)  as  well  as  the  path  of  integration  from  ooe~jn  to  oo.  In  general  there  also  exist 
a  finite  number  of  pole  singularities  (i.e.  N  number  of  poles)  on  the  top  Riemann  sheet; 
these  poles  are  called  proper  surface  wave  poles  (whose  residues  are  modal  waves  because 
they  satisfy  the  radiation  condition  at  infinity).  There  can  also  exist  an  infinite  number 
of  pole  singularities  on  the  bottom  Riemann  sheet;  the  latter  poles  axe  called  leaky  wave 
poles  (giving  rise  to  non-modal  waves  since  they  do  not  satisfy  the  radiation  condition  at 
infinity).  A  finite  number  of  surface  wave  poles  have  also  been  depicted  in  Figure  3.2. 

Equation  (3.34)  and  (3.35)  constitute  a  Sommerfeld  or  z-propagation  integral  repre¬ 
sentation;  on  the  other  hand,  the  p- propagation  representation  for  g  can  be  derived  by 
enclosing  all  the  singularities  of  g  in  the  plane  via  a  contour  deformation.  The  original 
contour  is  thus  deformed  as  shown  in  Figure  3.2  to  arrive  at 


u 2 

9i  r,r')  =  |2. 

4tt 


N 


(  g(Z,z,z')  nW(Sk0p)tdZ-2xj'£  Resmp,z,z'))ipH(02\t;pp) 

JCn 

(3.36) 


p=  i 


where  the  finite  sum  refers  to  the  discrete  spectral  or  pole  wave  contribution  to  g  and  the 
integral  of  Cg  is  the  continuous  spectral  contribution  to  g;  also,  the  Res  is  defined  as  the 
residue  of  g  at  the  poles  £  =  £p;  i.e.,  (see  also  (C.16)  and  (C.38)) 
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(3.37) 


Res{g(Zp,  z,  z '))  =  lim  (£  -  £p)  $((,  z,  z') 

s- %p 

The  integral  contour  Cg  can  be  transformed  to  the  real  axis  by  introducing  the  change  of 
variable 


;  <K  =  -jgs-  (3.38) 

It  is  noted  that  as  £  varies  along  the  contour  Cg  of  the  Figure  3.2,  the  new  variable  of 
integration  £  will  vary  from  +oo  to  — oo  along  the  real  £  axis  as  shown  in  Figure  3.3.  Note 
that  the  circular  contour  which  encircles  the  branch  point  does  not  contribute.  Finally,  one 
obtains 


g(T,  r')  =  H  $(£,*,*')  42)  (kopy/l-o)  £  d  £  (3.39) 

-2 *j  (^2Res(g((p,z,z)  £p  H^\^pk0p) 

The  above  result  in  (3.39)  represents  the  alternative  ^-propagation  representation  of  g  which 
is  obtained  from  its  Sommerfeld  representation  in  (3.34). 

The  radial  propagation  of  a  grounded  uniaxial  anisotropic  double  layered  dyadic  Green’s 
function  can  be  constructed  after  incorporating  g  of  (3.39),  U,  W,  and  V  given  in  the 
following  into  (3.14)  -  (3.17) 


(3.40) 


W<(r,  r')  = 


r  *§£$£<( C, *')  42)(ko pyjl  -  £2)£d£ 

J — oo  ^2 

-2*  j  £  ^rRes  \9'^p",z,z'))  H{02\koptp» 

p"  K2  \  ) 

t  {£  ^’(‘o,yrr 

77^~'Res  tfff'WV)) 

-  £  -77r-^es  Vv *1  *'))  #<S2)(*o  ?£P» 

n "  *2?P"  \  / 


(3.41) 
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(3.42) 


>  k3 

V<(r,rO  =  | 


*,  z1)  42\k0  pyfiTex  d  c 

J — oo  ^2 

~2njYl~rRes  (V«(£P',z,-0^  H^2)(k0pZP ') 


The  prime  (/)  and  double  prime  (//)  in  (3.40)  -  (3.42)  refer  to  TM  (or  E-mode;  a  -  1) 
and  TE  (or  H-mode;  a  =  2)  respectively;  see  also  Appendix  C.  One  can  now  see  that  the 
radial  propagation  representation  of  the  dyadic  Green’s  function  exhibits  fast  convergence 
for  laterally  separated  source  and  field  points.  This  is  due  to  the  fact  that  for  the  (  >  1, 
the  Hankel  functions  associated  with  the  integrand  in  (3.39)  -  (3.42)  decay  exponentially. 
Therefore  it  is  more  efficient  to  utilize  this  representation  for  the  planar  geometries  when 
dealing  with  a  solution  for  currents  within  a  planar  uniaxial  anisotropic  materials.  The 
logarithmic  singularity  of  the  Hankel  function  at  (  =  ±1  can  be  removed  via  the  following 
transformations 


<  =  y/l-t2  ,  ICI  <  1  ;  C  =  v/l  +  £2  >  ICI  >  1  •  (3.43) 

Hence,  the  final  form  of  g,  U,  W,  and  V  after  transforming  the  integrtal  (0,-1)  to  (0,1),  and 
(-oo,  -1)  and  (l,oo)  to  (0,  oo),  which  will  be  used  in  the  numerical  integration  are  given 
by: 


g(T,  r')  =  ^  [jf  [$(C,  -2,  A  -  A]  \(=v—Hq2)  (*o  pO  Z  d  Z 

+  [  _ _ Ho2)(-jkopO  Zd£ 

Jo  'c=v/l+£J 


~2*i  (e 

\  p 


($(£p,  z,  *)  Zp  rt2)(Spk0p) 


(3.44) 


*>,,')  =  h 


g't'fiC,  z,  A  -  ~9'u(-C,,z,  zr) 


[x  P2,t 

Jo  K'i(C) 

f°°  P2,t 

Jo  *'2'(C) 

-2 nj^^frRes  (g'f(tP»,z,z'))  H(2)(k0pZp„ 

p"  K2  \  / 


g£((;,z,z')-g'f(-<;,z,z') 


H^(k0pZ)ZdZ 

J|c=v/T^ 

42\-jk0pz)zdz 

j  lC=v/i+7T 


(3.45) 
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(0  ^d(Cz,z')-g'<(-Cz,z') 


P2  ,t 

*?(C) 


H[Q2)(k0pOj 

^c=Vi-«2 


-  f 


\  f2,t 


P2  ,t 

«?(C) 


HoH-jkopO^- 

^C=V,l  +  {2 


-2»>  ^E  7^1Res  ^P'’ z’  ^  #62'(fcopfp'))  ( 

-  E  4^:^"  (^(^p"’2’ z,))  ^o2)^°  ^p")j  |  * 


V<(r,r') 


1.3 

*o 

2tt 


1  e  2 


/  J 

JO 


(0 


E^XkopDldi 


c=>/i-<4 


too  £2 

>  > 

+  Jo  W 

-2niYl~rRes  ^.(V****^  Ho](koptp') 


H™(-jk0pOZdZ 

(3.47) 


The  behavior  of  the  integrand  of  the  radially  propagating  integral  representation  in  (3.44) 
to  (3.47)  will  be  discussed  in  Chapter  5  of  numerical  results. 

As  mentioned  earlier,  the  formal  representation  of  the  radially  propagating  integral  rep¬ 
resentation  of  the  dyadic  Green’s  function  developed  in  this  chapter  leads  to  a  very  efficient 
closed  form  asymptotic  approximation  of  the  Green’s  dyadic  which  remain  accurate  even 
for  very  small  lateral  separation  of  the  source  and  observation  points.  The  developement 
of  this  asymptotic  representation  will  be  discussed  in  the  next  chapter. 
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Figure  3.1:  Electric  point  current  dipole  source  in  a  grounded  anisotropic 
uniaxial  double  layered  slab,  (  e,  =  e,>tIt+eltizz  ;  pt,  =  /x,itIt+pifiZzz). 
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c‘ 


Figure  3.2:  Analytic  properties  of  g(£,z,z')  in  the  complex  £  plane;  sur¬ 
face  wave  poles  and  branch  cuts.  Also  depicted  here  are,  the  original 
Sommerfeld  contour,  the  contour  deformation  for  /^-representation  as  a 
sum  of  the  enclosed  residues  plus  the  integral  around  the  branch  cut 
Cb  of  gi^z^). 
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5" 


Figure  3.3:  Complex  C  plane  and  the  real  axis  contour  of  integration  of  (3.39). 
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Chapter  4 


Uniform  asymptotic  closed-form  representation  for  the  dyadic  Green’s 
function 


This  chapter  deals  with  the  development  of  an  asymptotic  representation  for  the  dyadic 
Green’s  function  for  the  uniaxial  anisotropic  double  layered  geometry  in  which  the  large 
parameter  directly  depends  on  the  lateral  separation  between  the  source  and  observation 
points;  however,  this  asymptotic  result  is  found  to  remain  accurate  even  for  source  and  field 
points  separated  by  only  a  few  tenths  of  a  free  space  wavelength.  As  mentioned  earlier, 
the  radially  propagating  integral  representation  given  in  the  previous  chapter  can  be  used 
as  a  starting  point  to  find  an  approximate  closed  form  asymptotic  solution  of  the  dyadic 
Green’s  function  which  remains  valid  in  the  entire  range  of  kQp  except  for  very  small  k0p. 
The  details  of  this  asymptotic  development  are  shown  in  the  following. 

Consider  the  integral  I  of  the  form 

1=  ^  F(C,z,z')H^\k0pJl-C2)  C  dC  ,  (4.1) 

in  which  ( kop )  is  assumed  to  be  a  large  parameter.  This  integral  is  similar  to  the  the  radially 
propagating  part  of  g,  U,  W,  and  V;  F(£,  z,  z')  has  the  following  functional  form 

F(<:,z,z')^K(Qe-jkoS<  ;  6  -  ±(z  ±  z')  (4.2) 

If  one  performs  the  angular  spectrum  mapping 

C  =  cos7;  dC,  = -sinydy,  and  \/ 1  -  C2  =  sin 7,  (4.3) 


then, 


e-jk0r  cos(0-7)  CQS  7  sin  7  dl  (4.4) 


where, 
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S  =  r  cos  9  ;  p  =  r  sin  9,  and  tan  9  =  |  , 


(4.5) 


and  use  has  been  made  of  the  large  argument  approximation  of  the  Hankel  function 

+  •  X>1  •  ,4-6) 


H(q2\x)  ~ 


The  integrand  exhibits  a  saddle  point  at  7,  =  6.  Let 


q(  7)  =  -jcos(0-7). 


(4.7) 


Then, 


q'il)  =  3  sin(0  -  7)  ;  7,  =  9  (saddle  point)  ;  q( jt)  =  -j  .  (4.8) 

Figures  4.1  and  4.2  show  the  double  sheeted  (  and  complex  7  planes  respectively,  with 
the  original  path  T  and  steepest  descent  path  Tsdp  depicted  in  the  complex  7  plane.  Note 
that  the  steepest  descent  path  is  the  path  on  which  the  following  condition  holds, 


!™{q{ 7))  =  Im(q( 7,)),  or  Re(cos{6  -  7))  =  1,  (4.9) 


so, 


f  f  1  +  ^  3k°rcos(0  1)  cos 7  sin 7^7 

Jr sdp  y  7rfc0rsin7  \  8  kQr  sin  7  / 


+  2 xj  Y,  Res(I( 7,))  ~  2x j  ^  Res(I(n)) 


(4.10) 


where  Tsdp  is  the  SDP  contour,  and  the  summation  terms  of  s  and  /,  stand  for  the  residue 
contributions  arising  from  the  presence  of  the  surface  wave  and  leaky  wave  poles  respectively 
which  may  be  intercepted  during  the  path  deformation.  It  is  noted  that  unless  8  «  -  (hop  » 
1;  k06  <  1),  there  exists  the  possibility  of  crossing  one  or  more  of  those  surface  wave  poles 
which  are  closest  to  the  Real  7  axis  (see  Figure  4.2)  which  have  been  captured  during 
the  derivation  of  p-propagation  representation  (see  (3.39)  -  (3.42)).  This  results  in  the 
cancellation  of  some  fo  the  surface  wave  residue  contributions  in  the  sum  on  the  right  hand 
sides  of  (3.39)  -  (3.42).  It  is  also  noted  that  the  same  SDP  will  result  if  one  starts  from  the 
Sommerfeld  representation  (rather  than  the  p-propagation  representation)  since  only  one 
SDP  exists  in  each  strip  of  27t  radians  in  the  7  plane,  but  in  this  case,  the  surface  wave 
poles  must  be  found  in  7  plane.  In  the  following  we  will  map  the  SDP  onto  a  contour  along 
the  real  axis  in  the  new  complex  plane  s  defined  by 
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9(7)- 9(7.)  =  s2  ;  cos(0  —  7)  =  1  —  js2, 


(4.11) 


or 


5  (branch  points  at  sjs2  +  2j  =  0).  (4.12) 

This  transformation  maps  the  steepest  descent  path  Tsop,  on  the  real  axis  of  the  complex 
3  plane  as  indicated  above.  The  saddle  point  at  7,  =  9  is  now  mapped  to  s  =  0  (see 
Figure  4.3).  For  the  6  ss  £  the  source  and  the  observation  points  laterally  separated  (i.e. 
( kop  1;  kQ6  «C  1),  one  will  have 


sin7  =  l-js2,  cos  7  =  |s|v/s2  +  2  j  . 


(4.13) 


It  is  useful  to  note  that 


C  =  sy/s2  +  2j  ;  f  = 


dC  2 j(s2  -  1) 


da  y/s2  +  2  j 

therefore  the  SDP  representation  of  I  in  the  s- plane  can  be  expressed  as 


(4.14) 


J~  G(s)e~ko'"7  ds-2njY  Res(I(i,))  ,  (4.15) 

y  ttkqP  J —00  i 

where,  in  view  of  (4.1),  G(s)  axe  given  by 

GW  =  0  +  *•  (416) 

It  should  be  noted  that  for  the  particular  case  of  6  w  ^  no  surface  wave  pole  will  be 
intercepted  during  the  path  deformation. 

The  uniform  asymptotic  approximation  of  J  can  be  represented  as,  [i9] 


7  e-]k°p  (  G"( 0) 


y/2je-—-  [ 

k0p  \  4  k0p 


-2*7  ^2  Rea(I(si))  , 

l 


(4.17) 
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where  s^'s  are  the  poles  (surface  wave  as  well  as  leaky  wave  poles)  that  are  considered 
to  be  close  to  the  saddle  point  at  s  =  0;  and  Res(G(s<,))  is  related  to  Res(F((<r))  (of 
p-representation)  and  Res(F(£„))  (of  the  Sommerfeld  representation)  by 


(4.18) 


and  s„  =  ±v/7\/^<r  -  1  !  C<r  =  ±\/l  -  also  s*  =  i/j/V1  ~Q~  1-  The  value  of 
G"(s),_0  is  related  to  the  first  derivative  of  F(s)  by 


G»'-°  ■  «‘  +  Sih>Z'Wi-  . 


=  _4j<1  +  siw)v^7  ’ 


(4.19) 


where  use  has  been  made  of  gj  =  and  £  =  1  —  js 2  ;  =  —2 js.  The  transition 

function  .F(a)  is  defined  [19,  20] 


F(a)  -  2 jyjaeic*  f 


e~>r  dr  , 


(4.20) 


the  analytical  property  of  the  transition  function  is  also  discussed  in  [19]. 
The  function  F(£)  can  in  general  be  represented  as 

evc  -  J\  _  zi  Z  ) 


=  m- 


D(0  ’ 


(4.21) 


where 


v/w2  /(£)]  =  0  ;  lim  [v/l  -  £2  ^/(o]  =  0  .  (4.22) 


lira  yfi ±-F(Z,z,z') 


m  N(Z,  z,  z')  +  f(i)  N'{S,z,z') 


=  limv/TT^ 

?  15(0 


i?(l)  U 


Z?(l)  « 


(4.23) 


lim  v/l  -  £2  D'((),tml 
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Hence  the  uniform  asymptotic  representation  of  a  grounded  uniaxial  anisotropic  double 
layered  dyadic  Green’s  function  can  be  constructed  after  incorporating  (4.17)  (in  view  of 
(4.18),  (4.19),  (3.39)  -  (3.42))  into  (3.14)  -  (3.17),  where  fy(£),fy(0.* yytf),  and  Fv(0 
used  in  (4.18)  and  (4.19)  are  given  in  the  following: 


Fg<{Z,z,z') 

FU(t,z,z') 

Fw(S,z,z’) 


De(0  ’ 


=  fci 


^2 ,t  Z'  2  ) 

De(  0  ’ 

k0  (  k2  z,  z')  n2  t  vIh(Z,z,z')) 


?  *2 ,t  DE(0 


Dh(  0 


(4.24) 


(4.25) 


(4.26) 


Fv(t,z,z’)  =  k^ 


e  N<<e(^z,z') 


*2  De(0 


(4.27) 


and  expressions  for  N$( l),limj 


H 


,  DE(l),lim£_1 
H 


Vf-T*  1£De(0 


H 


and  the  residues  of  g[f(z,  z')  and  g^(z,z')  are  all  given  in  Appendix  C. 


The  uniform  asymptotic  closed  form  approxiamtion  of  the  dyadic  Green’s  function  de¬ 
veloped  here  has  been  cast  in  such  a  form  that  the  physical  behavior  of  the  electromagnetic 
fields  due  to  anisotropy  of  the  medium  reveals  itself  through  a  simple  mathematical  pa¬ 
rameters  (see  also  Appendic  C).  Moreover,  the  asymptotic  solution  of  each  of  the  g( r,r'), 
W(r,  r'),  and  V(r,  r')  consists  of  the  contribution  of  the  space  wave  (the  saddle 
ponit  contribution),  surface  waves,  leaky  waves  (if  any),  and  the  transitoin  effects  of  the 
surface  and  leaky  wave  poles  which  axe  close  to  the  saddle  point.  Thus,  the  physical  under¬ 
standing  of  the  interaction  of  the  spatially  confined  source  with  an  anisotropic  (uniaxial) 
double-layered  grounded  slab  is  greatly  enhanced  through  the  newly  developed  asymptotic 
closed  form  representation  of  the  dyadic  Green’s  function. 


Becuase  there  are  infinite  number  of  leaky  wave  poles  in  the  bottom  sheet  of  the  Som- 
merfeld  plane;  the  question  may  arise  which  poles  should  be  cosiderd  close  to  the  saddle 
point.  For  the  loss-less  case,  there  are  two  sets  of  leaky  wave  poles  on  the  bottom  sheet 
of  the  Sommerfeld  plane;  the  complex  poles  and  the  real  poles.  During  the  investigation  it 
is  found  that  the  inclusion  of  the  real  poles1,  which  may  be  called  improper  surface  wave 
poles  [20,  23,  24],  (because  they  are  the  same  species  as  of  surface  wave  poles  lying  on  the 

'They  are  lying  on  the  positive  imaginary  axis  of  {  plane,  the  radially  propagating  plane. 
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bottom  sheet  of  the  Sommerfeld  plane),  would  be  sufficient  for  achieving  a  good  accuracy 
for  the  asymptotic  representation. 

In  fact  the  proper  surface  wave  poles  are  originally  improper  surface  waves  [29,  24]: 
of  course  except  for  the  first  TM  (E-mode)  surface  wave  mode  which  is  always  present 
independent  of  material  properties  and  thicknesses  of  the  layers.  As  the  materials  become 
more  dense  and/or  the  substrates  become  electrically  thicker,  an  improper  surface  wave 
pole  moves  toward  the  brach  point,  and  it  eventually  passes  through  the  branch  point  and 
appears  in  the  upper  sheet  of  the  Sommerfeld  plane;  and  hence,  it  becomes  a  proper  surface 
wave  pole.  It  is  easy  to  verify  that  the  condition  for  any  surface  wave  mode  to  “turn-on”  is 
£p  =  1.0,  (note  that  the  variable  f  has  been  normalized  with  respect  to  k0).  From  (C.12) 
for  TM  (E-mode),  one  will  have 


-De(£p)|(p,i  0  -  0  , 


(4.28) 


or 


A-  tan  L\  = - A  tan  L'2  ,  (4.29) 

f2,t“2 

where  L\  =  diy/n,  -  n' ,  n,  =  and  d,  =  k0  d,.  By  the  same  token,  from 

(C.34)  for  TE  (H-mode),  one  will  have 

c»(f,)|4., ,  =  0  ■  («<>> 


or 


L'i 


tan  L'[ 


L" 

^2,^2 


cot  L'2 


(4.31) 


where  L"  =  d,  yj n,  -  n'\  and  n "  = 

Simillar  equations  for  the  isotropic  case  have  already  been  reported  by  Alexopoulos 
and  Jackson  [25].  Equations  (4.29)  and  (4.31)  respectively  represent  parametric  conditions 
for  any  T M  and  T E  surface  wave  modes  to  turn-on.  The  loci  of  these  equations  for  the 
isotropic  case  have  been  studied  by  Marin  and  Pathak  [26].  It  is  evident  that  by  controlling 
different  parameters  invloved  in  (4.29)  and  (4.31),  one  can  excite  or  suppress  a  certain 
number  of  poles.  A  complete  elimination  of  the  surface  wave  excitation  is  also  possible  if 
one  lets  k'2(£p)  =  0  for  TM  (E-mode),  and  ^((p)  =  ®  ^or  TE  (H-mode).  This  is  due  to 
the  fact  that  since  there  is  no  variation  of  the  surface  wave  modes  in  2-direction  with  these 
conditions,  and  the  tangential  electric  field  must  vanish  on  a  perfectly  conducting  ground 
plane;  by  reciprocity  this  impies  that  these  modes  are  not  excited  [25]. 
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It  is  also  important  to  realize  that  because  the  proper  and  improper  surface  wave  poles 
are  lying  on  the  real  axis  of  the  Sommerfeld  plane  (positive  imaginary  axis  of  the  radially 
propagating  plane),  they  can  easily  be  found2;  and  hence,  rendering  to  the  efficeincy  of  the 
newly  developed  asymptotic  solution. 

In  the  following  chapter  the  accuracy  and  the  efficiency  of  the  asymptotic  solution  as 
a  function  of  lateral  separation  of  the  soure  and  observation  points  with  different  material 
parameters  and  thicknesses. 


2  A  Newton- Raphson  search  procedure  with  a  good  initial  guess  can  be  employed  to  find  the  poles  in  a 
few  iterations. 
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Figure  4.1:  Two  sheeted  £  plane.  First  sheet,  —  £2)  <  0;  and  the 

second  sheet,  Im(y/ 1  -  £2)  >  0. 


48 


LEAKY  WAVE  S.W.  POLES 


r,„  r  poles 


Figure  4.2:  Angular  spectrum  mapping  where  (  =  cos 7;  the  leaky  and 
surface  wave  poles  and  the  contour  of  integration  T  and  Tsdp 
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Figure  4.3:  Complex  s  plane,  original  path  T  maps  on  solid  curve,  the  SDP 


maps  onto  the  real  axis  (—00  <  s  <  -00),  the  saddle  point  7,  =  0,  and 
branch  points  at  =  ±y/2 e~^ir^4(^s2  +  2 j  =  0). 
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Chapter  5 


Numerical  Results 


Figure  5.1  shows  the  typical  behavior  of  the  integrand  of  the  function  W  of  (3.31)  for  p/  A0  = 
0.1  versus  £  as  \z  -  z'\  — >  0.  The  integrand  is  an  oscillatory  function  (typical  behavior  of  the 
Bessel  function),  and  decays  very  slowly.  Figure  5.2  shows  the  numerical  behavior  of  the 
integrand  of  W  of  (3.33)  for  p/A0  =  0.1  as  \z  -  z'\  — *  0  when  the  singularity  of  the  integrand 
is  removed.  The  significance  of  the  singularity  removal  technique  is  vivid  if  one  compares 
the  smoothness  of  the  plots  in  Figures  5.1  and  5.2  for  p/ Ao  —  0.1.  It  is  now  evident  that  (for 
the  small  lateral  separation)  the  method  of  singularity  removal  technique  can  efficiently  be 
utilized  to  evaluate  the  dyadic  Green’s  function.  Nevertheless  as  the  lateral  separation  of  the 
source  and  observation  points  increases  (kop  becomes  large)  the  integrand  starts  to  oscillate. 
Figure  5.3  shows  the  typical  behavior  of  the  integrand  of  W  of  (3-33)  for  pj  Ao  =  4.7.  This 
oscillation  is  due  to  rapid  growth  of  the  argument  of  the  Bessel  function  which  causes  the 
method  to  lose  its  efficiency.  This  is  the  main  difficulty  of  the  singularity  removal  technique 
which  motivates  us  to  construct  the  alternative  radial  propagation  representation  of  the 
dyadic  Green’s  function.  Figures  5.4  and  5.5  show  a  typical  behavior  of  the  integrand 
of  W  of  (3.46)  as  a  function  of  its  argument  £  for  different  lateral  distances  p/A0  =0.5, 
and  2.0.  It  is  important  to  note  that  as  the  separation  between  the  source  and  observation 
points  increases,  the  infinite  integral  decays  faster,  but  the  number  of  oscillations  in  the 
finite  integrals  increases.  This  can  be  seen  from  Figure  5.5.  These  oscillation  are  due  to 
the  behavior  of  the  Hankel  function  as  its  argument  becomes  larger.  It  is  interesting  to 
note  that  in  the  Sommerfeld  representation,  the  total  number  of  oscillation  of  the  integrand 
between  0  and  1  of  its  semi-infinite  interval  of  integration  are  the  same  as  the  number  of 
oscillations  of  the  integrand  of  the  radial  propagation  representation  which  exist  only  over 
the  finite  limits  from  0  to  1.0.  Over  the  rest  of  the  interval  of  the  integration  from  1.0 
to  infinity,  however,  the  Sommerfeld  integrand  exhibits  an  infinite  number  of  oscillations; 
no  such  oscillations  exist  for  the  radial  propagation  representation.  The  advantage  of  the 
radial  propagation  representation  against  the  Sommerfeld  representation  is  now  apparent. 

Figures  5.6-5.13  show  the  numerical  integration  and  the  closed  form  asymptotic  repre¬ 
sentation  of  W  given  in  (4.17)  versus  p/A0  for  various  substrate  and  superstrates  thicknesses 
with  high  and  low  material  constants  and  various  degree  of  anisotropy.  The  numerical  val¬ 
ues  of  these  parameters  are  given  in  the  captions  of  the  Figures.  The  accuracy  of  the  new 
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representation  is  quite  impressive.  This  new  asymptotic  solution  remains  valid  for  the  field 
points  very  close  to  the  source.  The  efficiency  and  hence  its  usefulness  of  the  new  asymp¬ 
totic  expansion  is  now  clear.  The  error  bound  of  the  solution  is  less  than  one  to  five  percent 
for  lateral  separations  between  the  source  and  field  points  as  small  as  0.5- 1.0  of  free  space 
wavelength.  Also  it  is  note  worthy  that  as  the  lateral  separation  of  the  source  and  obser¬ 
vation  increases,  the  computational  time  for  the  numerical  integration  also  increases.  This 
time  increase  is  due  to  rapid  growth  of  the  argument  of  the  Bessel  and  Hankel  functions 
which  are  present  in  the  integrand.  Nevertheless,  the  computational  time  for  the  asymp¬ 
totic  representation  remains  constant  since  no  Sommerfeld  type  integral  is  involved  in  this 
solution;  and  moreover,  the  accuracy  of  the  asymptotic  solution  improves  as  the  lateral 
separation  of  the  source  and  observation  point  increases. 

As  expected,  it  is  found  that  the  accuracy  of  the  new  representation,  particularly  for 
small  lateral  separation  of  the  source  and  observation  points,  indeed  depends  on  the  mate¬ 
rial  thicknesses  and  parameters.  In  the  case  of  thin  slabs  with  low  material  parameters,  the 
dominant  effect  is  due  to  space  waves  arising  from  the  saddle  point  contribution,  whereas 
for  thick  and  dense  material  constants,  the  dominant  effect  is  due  to  the  surface  waves.  For 
the  uniform  asymptotic  approximation  only  one  TM  surface  wave  pole  to  be  considered 
close  to  ‘He  saddle  point  for  thin  and  low-dense  material  slabs.  For  thicker  and/or  more 
dense  material  slabs,  a  TE  improper  surface  wave  pole  also  comes  close  to  saddle  point, 
and  as  discussed  in  the  previous  chapter,  when  a  certain  condition  met,  (see  (4.31)),  the 
TE  improper  surface  pole  eventually  becomes  the  first  TE  proper  surface  wave.  In  most 
practical  applications,  the  uniform  asymptotic  dyadic  Green’s  function  needs  one  TM  sur¬ 
face  wave  and  one  T  F  proper  or  improper  surface  wave  pole  singularities  close  to  the  saddle 
point.  It  is  noted  ‘hat  the  improper  surface  poles  are  located  on  the  positive  imaginary 
axis  of  the  £  plane  ^ radially  propagating  plane),  which  are  the  same  improper  surface  wave 
poles  lying  on  the  positive  real  axis  of  the  bottom  sheet  of  the  original  Sommerfeld  plane. 

The  accuracy  of  the  uniform  asymptotic  representations  of  g( r,  r'),  U{ r,  r'),  and  V(r,  r') 
have  been  tested,  and  they  are  at  least  as  accurate  as  of  W(r,  r').  Those  results  have  not 
been  reported  here. 

A  comparison  between  the  numerical  and  asymptotic  evaluation  of  electric  fields,  Ex,  Ey, 
of  an  electric  point  current  parallel  to  the  interface  (horizontally  oriented)  over  the  dielectric 
interface  for  the  grounded  double- layered  anisotropic  uniaxial  slab  versus  <f>  at  p  =  0.5Ao 
are  shown  in  Figures  5.14  -  5.17.  It  can  be  seen  from  the  figures  that  the  accuracy  of  the 
asymptotic  representation  is  quite  accurate  and  remains  valid  even  for  the  field  points  very 
close  to  the  source. 

With  the  new  asymptotic  solution  the  effect  of  the  material  anisotropy  over  the  field 
components  can  easily  be  investigated.  The  parameters  n'  =  and  n"  =  ^  represent 
the  degree  of  anisotropy  of  the  material,  and  can  have  any  value  in  complex  plane.  It 
should  be  noted  that  as  n'  and  n"  go  to  one  the  m^*erial  becomes  isotropic.  Also  if  the 
thickness  of  the  first  layer  goes  to  zero,  or  its  material  parameters  approach  to  the  free 
space  constitutive  parameters  the  solution  will  give  the  dyadic  Green’s  function  Tv  a  single 
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layer  grounded  anisotropic  uniaxial  slab. 

The  effect  of  the  anisotropy  on  the  W(r,  r')  ,  can  be  seen  if  one  compares  the  plots  of 
Figure  5.6  (for  the  isotropic  case),  and  Figure  5.7  (for  n[  —  n'2  =  1.5),  and  Figure  5.8 
(for  n'x  =  n2  =  0.75).  It  seems  that  the  strength  of  the  field  increases  as  the  value  of  nj 
decreases.  Similar  conclusion  can  be  drawn  from  the  plots  of  the  surface  electric  fields  of 
shown  in  Figures  5.14  -  5.16. 

More  numerical  results  and  the  physical  interpretations  based  on  the  newly  developed 
closed  form  asymptotic  representation  of  the  dyadic  Green’s  function  for  a  double-layered 
anisotropic  uniaxial  grounded  slab  will  be  reported  in  the  second  phase  of  this  research. 
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£  =  ALPHA 


p  =  0.1  Xq 


Figure  5.1:  Behavior  of  the  integrand  of  the  Sommerfeld  type  integral  of 
W  (3.31)  for  p/Xo  =  0.1  versus  £  as  \z  -  z'\  — ►  0. 
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Figure  5.2:  Behavior  of  the  integrand  of  the  Sommerfeld  type  integral  of 
W  in  (3.33)  for  p/X0  =  0.1  versus  f  as  \z-z'\  -*•  0,  when  the  singularity 
removel  technique  is  employed 
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P  =  4.7  A0 


Figure  5.3:  Behavior  of  the  integrand  of  the  Sommerfeld  type  integral  of 
W  in  (3.33)  for  p/XQ  =  4.7  versus  £  as  \z-z'\  — *•  0,  when  the  singularity 
removel  technique  is  employed 
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p- representation 

- Real  part 

- Imaginary  parts 


p  =  0.5Ao 


Figure  5.4:  Behavior  of  the  integrand  of  the  radially  propagating  type 
integral  of  W  in  (3.46)  for  p/A0  =  0.5  versus  £  as  \z  -  z']  —  0. 


57 


p- representation 

- Real  part 

- Imaginary  parts 


p  =  2.A0 


Figure  5.5:  Behavior  of  the  integrand  of  the  radially  propagating  type 
integral  of  W  in  (3.46)  for  p/XQ  =  2.0  versus  £  as  \z  -  z'\  -*  0. 
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W  (Real  &  Imaginay) 


•05,  eM  -  2.25,  clit/c1  z  =  1.0,/ii,t  =  1.0,  fti,t/fii>z  =  1.0 
•05,  e2,t  =  3.25,c2,j/e2  ,  =  1.0,  ^2,t  =  1.0,/r2it/^2vt  =  1.0 
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Figure  5.6:  Comparison  between  the  numerical  integration  and  asymp 
totic  value  of  W  versus  p/A0  in  (4.17)  for  |z  -  z'\  0. 


d\/ \q  —  0.05,  ei.t  =  2.25,  eu/ei,*  —  l-5»  A^i,*  —  1-0,  Pi,t/ f*i,x  1-0 
J2/A0  =  0.05,  t2,t  —  3.25,  62,t/e2 ,1  =  1-5, /^2,t  =  1-0,  Ai2,t/^2,z  =  1-0 


Figure  5.7:  Comparison  between  the  numerical  integration  and  asymp¬ 
totic  value  of  W  versus  pj Aq  in  (4.17)  for  \z  —  z'\  — ►  0. 
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W  (Real  &  Imaginay) 


d\/\o  —  0.05,  ex,*  —  2.25,  Ci.t/ej,,  =  0.75,  =  1.0,  Pi,t/ Pi,*  —  1-0 

dz/Xo  =  0.05,  €2, t  =  3.25,  e2,t/€2,i  =  0.75,  ^2, t  =  1-0,  A*2,t/M2,*  =  1-0 


(p/x 0) 


Figure  5.8:  Comparison  between  the  numerical  integration  and  asymp¬ 
totic  value  of  W  versus  p/Xq  in  (4.17)  for  \z  —  z'\  — ►  0. 
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W  (Real  &  Imaginay) 


5,Ci,t  —  3.25,  ei.i/ei,,  =  1.0,  H\,t  =  1.0,  =  10 

3,€2it  =  9.60,e2,4/e2,i  =  1.0,  /i2, 4  =  1.0,/42>4/^2,x  =  1.0 
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.  Im.  (Asymptotic) 

Figure  5.9:  Comparison  between  the  numerical  integration  and  asymp 
totic  value  of  W  versus  p/X 0  in  (4.17)  for  \z  —  z'\  — »  0. 


W  (Real  &  Imaginay) 


dxf  A0  =  0.05,  €i,e  =  3.25,  =  1.5,/n.*  =  1.0,/W/*M  =  10 

d2/X o  =  0.03,  e2,t  =  9.60,  e2,t/e2ii  =  0.9,  M2, t  =  1-0,  /W M2.*  =  1-0 
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Figure  5.10:  Comparison  between  the  numerical  integration  and  asymp¬ 
totic  value  of  W  versus  p/ Aq  in  (4.17)  for  \z  —  z’\  — ♦  0. 
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W  (Real  &  Imaginay) 


dJXQ  =  0.05, 6i,t  =  3.25,«w/eM  =  =  WPulPu  ~  10 

d2/X0  =  0.03,  £2.t  =  9.60  ,e2it/e2>z  =  0.5,/i2,t  =  1.0,fi2,t/p2,z  =  10 
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Figure  5.11:  Comparison  between  the  numerical  integration  and  asymp¬ 
totic  value  of  W  versus  p/X o  in  (4.17)  for  \z  —  z'\  — ►  0- 
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W  (Real  &  Imaginay) 


di/ Ao  =  0.05,  ei>t  =  3.25,  fi,*/ tj,,  —  1.0, mm  —  1-0, Mm/ Mi,*  —  1-0 
d2/ A0  =  0.03,  e2,i  =  2.25,  e2,t/e2,*  =  1-0,  M2,  t  =  1-0,  p2,t/ P2,x  =  1-0 
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Figure  5.12:  Comparison  between  the  numerical  integration  and  asymp¬ 
totic  value  of  W  versus  pf  Aq  in  (4.17)  for  \z  —  z'\  —*  0. 
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W  (Real  &  Imaginay) 


d\/  V)  —  0.03,  £iit  =  2.25,  £i,t/ £jiZ  =  1.0,  /iiit  =  1.0,  pi,t/p i,z  =  1.0 
^2/^0  =  0.06,  62, t  =  9-60,  t2,t/e2,i  =  10,  p2,t  =  1.0,  P2,t/P2,z  =  1-0 


Figure  5.13:  Comparison  between  the  numerical  integration  and  asymp 
totic  value  of  W  versus  p/Xo  in  (4.17)  for  | z  -  z'\  — ►  0. 
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di/\o  —  0.05,  €ift  —  2.25,  ciit/cit,  =  =  1-0,  pi,t/ Pi,z  =  1.0 

rf2/A0  =0.05,e2,t  =  3.25,  e2,t/e2ii  =  1.0,/x2if  =  1.0,  ^2,t//i2,z  =  1.0 


Figure  5.14:  Comparison  between  the  numerical  integration  and  asymp¬ 
totic  value  of  electric  fields  Ex  and  Ey  of  a  horizontal  electric  dipole  cur¬ 
rent  point  source  over  the  planar  interface  of  the  two  anisotropic  uniax¬ 
ial  slab  at  observation  points  (p,  4> )  versus  <f>  at  (p  =  0.5  Ao,  \z  —  z'\  =  0) 
after  incorporating  (4.17)  into  (3.14)-(3.17). 
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5,€i,t  —  2.25,  ei(t/e1>z  =  . 

L-5,  H\jt  =  1-0,  =  1.0 

5,  e 2,t  =  3.25,  62, t/ ^2,*  = 

L.5  =  1.0,p2,t/A‘21*  =  1-0 

■  Sommerfeld 

(p  =  0.5  A0,  | z  -  z'\  =  0) 

<  Asymptotic 

Figure  5.15:  Comparison  between  the  numerical  integration  and  asymp 
totic  value  of  electric  fields  Ex  and  Ev  of  a  horizontal  electric  dipole  cur 
rent  point  source  over  the  planar  interface  of  the  two  anisotropic  uniax 
ial  slab  at  observation  points  (p,  <£)  versus  <f>  at  (p  =  0.5  A0, \z  —  z'\  =  0 
after  incorporating  (4.17)  into  (3.14)-(3.17). 


<ii/Ao  =  0.05,  ei,t  =  2.25,  ,*/ ei,*  —  0.75,  pi>t  =  1.0,  pi ,t/pi,z  =  1-0 
^2/^0  =  0.05,  €2, t  =  3.25,  €2,t/c2,z  —  0.75, /*2 ,t  =  1-0,  P2,t/^2,z  =  1-0 
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Figure  5.18:  Comparison  between  the  numerical  integration  and  asymp¬ 
totic  value  of  electric  fields  .Ex  and  Ey  of  a  horizontal  electric  dipole  cur¬ 
rent  point  source  over  the  planar  interface  of  the  two  anisotropic  uniax¬ 
ial  slab  at  observation  points  (p,  <f> )  versus  4>  at  (p  =  0.5  Ao,  | z  —  z'\  =  0) 
after  incorporating  (4.17)  into  (3.14)-(3.17). 
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^iAo  =  0.05,  €i,t  =  3.25,  ei.t/ei,*  =  0.5,piit  =  1.0,  pi.t/pi,*  =  1.0 
^2/^0  =  0.03,  €2,t  =  9.60,  €2,t/€2,i  =  0.5,  fi2,t  =  1.0,  /i2,t/A£2,z  =  10 


Figure  5.17:  Comparison  between  the  numerical  integration  and  asymp¬ 
totic  value  of  electric  fields  Ex  and  Ey  of  a  horizontal  electric  dipole  cur¬ 
rent  point  source  over  the  planar  interface  of  the  two  anisotropic  uniax¬ 
ial  slab  at  observation  points  (p,  <£)  versus  <j>  at  (p  =  0.5  Ao,  \z  —  z'\  =  0) 
after  incorporating  (4.17)  into  (3.14)-(3.17). 
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Chapter  6 


Conclusions,  Discussions  and  Directions  of  Future  Research 


A  relatively  simple  and  accurate  closed  asymptotic  form  of  the  electric  dyadic  Green’s  func¬ 
tion  for  a  double-layered  anisotropic  uniaxial  grounded  slab  is  developed.  The  asymptotic 
representation  provides  an  accurate  solution  even  for  a  very  small  separation  of  the  source 
and  observation  points  (a  few  tenths  of  a  free  space  wavelength).  The  asymptotic  closed 
form  dyadic  Green’s  function  has  been  cast  in  such  a  form  that  the  physical  behavior  of 
the  electromagnetic  fields  due  to  anisotropy  of  the  medium  reveals  itself  through  a  simple 
mathematical  parameters.  Thus,  the  physical  understanding  of  the  interaction  of  the  spa¬ 
tially  confined  source  with  an  anisotropic  (uniaxial)  double-layered  grounded  slab  is  greatly 
enhanced  through  our  asymptotic  closed  form  representation  of  the  dyadic  Green’s  func¬ 
tion.  In  addition  two  efficient  numerical  integration  schemes  based  on  singularity  removal 
technique  and  radially  propagating  integral  representation  of  the  Green’s  dyadic  are  also 
presented.  In  contrast  to  the  representations  developed  in  this  study,  the  Sommerfeld  type 
integral  converges  slowly  for  large  and  even  moderately  large  separation  between  the  source 
and  observation  points.  The  representations  developed  here  do  not  exhibit  this  difficulty. 
Hence  a  systematic  numerical  algorithm  is  developed  which  can  efficiently  be  used  to  com¬ 
pute  the  dyadic  Green’s  function  of  the  anisotropic  (uniaxial)  double-layered  grounded  slab 
for  the  entire  range  of  the  separation  of  source  and  observation  points.  These  solutions  are 
expected  to  be  useful  in  the  moment  method  solution  to  a  wide  variety  of  the  anisotropic 
uniaxial  grounded  slab  problems. 

Some  related  solutions  such  as  the  dyadic  Green’s  function  for  electric  as  well  as  magnetic 
point  currents  in  a  general  uniaxial  multilayered  material  media  have  been  developed  in  a 
form  which  provide  useful  physical  interpretations;  these  solutions  are  important  in  many 
applications  involving  uniaxial  stratified  media,  such  as  in  the  design  of  high  performance 
large  phased  arrays  and  frequency  selective  surfaces  (FSS),  as  well  as  in  the  analysis  of 
multilayered  guided  wave  structures. 

It  is  believed  that  the  derivations  of  the  Green’s  function  for  multi-layered  anisotropic 
(uniaxial)  media,  as  well  as  the  efficient  asymptotic  closed  form  approximation  to  the  dyadic 
Green’s  function  for  anisotropic  (uniaxial)  double- layered  grounded  slab  are  novel;  and  be¬ 
cause  it  is  an  efficient  solution  to  a  canonical  problem,  it  has  applications  in  a  variety  of 
disciplines;  such  as,  monolithic  microwave  and  millimeter  waves  integrated  circuits,  inte- 
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grated  optics,  non-destructive  evaluation  (NDE),  and  geophysical  prospecting.  As  men¬ 
tioned  earlier,  the  conventional  spectral  integral  of  the  Green’s  dyadic  is  highly  oscillatory 
and  contains  branch  points  and  pole  singularities  that  are  on,  or  close  to,  the  path  of  in¬ 
tegration.  Hence,  the  numerical  evaluation  of  the  spectrum  integrals  is  in  general  very 
time  consuming  and  becomes  prohibitively  expensive  as  the  separation  of  the  source  and 
observation  points  increases  in  the  lateral  direction. 

The  utility  of  our  efficient  asymptotic  closed  form  approximation  for  the  Green’s  dyadic 
can  hardly  be  overemphasized. 

During  the  Phase  II  of  our  research  effort  we  will  perform  the  following  : 

a)  During  Phase  I  of  our  work,  we  have  gained  a  great  deal  of  physical  insight  into  the  elec¬ 
tromagnetic  source  interaction  with  anisotropic  media.  Owing  to  the  closed  form  expression 
for  the  asymptotic  dyadic  Green’s  function  for  anisotropic  double-layered  grounded  slab, 
the  effect  of  anisotropy  reveals  itself  through  different  mathematical  parameters.  Because 
of  the  presence  of  ten  parameters  (five  for  each  anisotropic  layer  including  the  thickness 
of  the  layer)  extensive  numerical  simulation  must  be  performed  to  completely  character¬ 
ize  a  source  excited  anisotropic  double-layer  geometry.  We  will  perform  such  an  extensive 
numerical  simulation  during  Phase  II. 

b)  We  will  extend  our  dyadic  Green’s  function  analysis  to  include  the  Greeu’s  dyadic  due  to 
magnetic  point  current.  This  analysis  allows  us  to  expand  the  applicability  of  our  method 
to  aperture  coupled  and  co-planar  integrated  circuit  devices. 

c)  We  intend  to  develop  a  computer  aided  engineering  (CAE)  code  to  analyze  monolithic 
microwave  and  millimeter  wave  integrated  planar  elements  as  shown  in  Figures  6.1,  6.2,  and 
6.3.  This  will  include:  1)  the  analyses  of  passive  microwave  and  millimeter  wave  integrated 
circuits  such  as  design  of  arbitrary  bends,  discontinuities,  filters,  and  couplers,  2)  finite 
as  well  as  infinite  microstrip  antenna  phased  array  (including  the  feeding  network),  3) 
integrated  optical  devices  such  as  dielectric  waveguide  with  arbitrary  bends,  and  dielectric 
waveguide  couplers  all  in  anisotropic  uniaxial  double-layered  grounded  slab  environment. 

d)  We  will  develop  a  computer  code  to  analyze  the  integrated  optical  components  which 
needs  a  calculation  of  the  volumetric  current  distribution.  To  be  more  specific,  we  will  ana¬ 
lyze  a  dielectric  waveguide,  with  arbitrary  bends,  residing  within  an  anisotropic  (uniaxial) 
double-layered  grounded  slab.  This  configuration  is  a  canonical  high  frequency  transmission 
channel  of  the  electromagnetic  energy,  and  its  analysis  is  of  great  interest  to  the  optical  com¬ 
munity.  It  is  noteworthy  that  because  we  have  an  efficient  Green’s  dyadic  of  the  anisotropic 
double-layered  structure,  the  only  unknown  we  have  is  the  volumetric  current  within  the 
transmission  channel,  which  we  will  solve  for  by  our  efficient  conjugate  gradient  -  fast 
Fourier  transform  (GC-FFT)  method  that  we  will  explain  shortly.  Nevertheless,  because 
of  the  presence  of  the  open  boundaries  the  finite  element  analysis  of  this  class  of  problems 
is  very  time  consuming;  it  takes  perhaps  (10-100)  times  more  computational  effort  than 
our  volume-integral  formulation  which  utilizes  the  newly  developed  efficient  dyadic  Green's 
function  as  a  kernel. 
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e)  In  the  past,  we  have  used  the  unit  pulse  function,  defined  on  a  regular  three-dimensional 
grid,  to  discretize  the  volume-integral  equation.  The  advantage  of  this  discretization  scheme 
is  that  it  creates  a  resulting  matrix  Toeplitz-Hankel,  which  greatly  facilitates  the  compu¬ 
tation  and  storage  demands  on  the  computer.  Pulse  functions  however  do  not  provide 
continuity  of  the  currents  from  one  cell  to  another,  so  it  is  desirable  to  consider  other  func¬ 
tions  that  do  provide  this  continuity.  Recently,  we  have  introduced  splines  that  are  derived 
from  higher-order  convolutions  of  the  unit  pulse,  which,  of  course,  are  also  defined  on  the 
same  regular  grid  as  are  the  pulse  functions.  We  have  also  introduced  the  reciprocal  bases 
over  the  same  regular  grid,  which  are  orthogonal  to  the  spline  functions.  When  we  use  our 
higher-order  splines  as  expansion  functions  and  the  associated  reciprocal  bases  as  a  testing 
function  to  discretize  the  volume-integral  equation,  we  will  retain  the  desirable  feature  of 
the  Toeplitz-Hankel  matrix  in  our  system  of  linear  equations  as  well  as  the  continuity  of 
the  currents. 

It  is  also  possible  to  use  higher-order  shape  functions,  that  are  typically  used  in  the 
finite-element  method,  to  discretize  the  volume-integral  equation.  The  advantage  of  higher- 
order  shape  functions  is  that  they  not  only  provide  the  continuity  of  sources,  but  also 
they  can  easily  conform  to  the  shape  of  the  scatterer,  and  thereby  give  a  more  accurate 
discretized  version  of  the  continuum.  However,  due  to  the  irregular  three-dimensional  grid 
producing  by  this  discretization  scheme,  the  resulting  matrix  may  not  be  in  Toeplitz  or 
Hankel  form. 

We  will  utilize  both  discretization  techniques  (higher-order  splines,  and  higher  order 
shape  functions)  for  our  volume-integral  equation.  We  will  study  the  advantages  of  each 
method  to  our  problem.  We  will  also  look  into  the  use  of  combination  of  both  discretization 
methods  to  take  advantage  of  both  discretization  methods.  This  can  be  done  by  using  the 
finite  element  shape  functions  only  in  the  boundary  regions  and  the  higher-order  splines 
in  the  rest  of  the  body  of  scatterer.  This  hybrid  method,  although  it  leads  to  a  complex 
discretization  procedure,  is  expected  to  produce  more  accurate  results  than  the  original 
higher-order  spline  technique  that  makes  an  staircase  approximation  of  the  outer  boundary 
of  the  scatterer. 

We  will  also  study  the  degree  of  the  complexity  of  the  discretization  procedure  and 
make  a  decision  in  the  early  stages  of  Phase  II  as  to  which  method  is  more  suitable  for  the 
discretization  of  our  volume- integral  equation. 

f)  Conjugate  gradient  method  and  its  pre-conditional  version  is  a  very  powerful  tool 
to  solve  a  system  of  matrix  equations.  When  it  can  be  combined  with  the  Fast  Fourier 
Transform  (FFT)  algorithm,  it  will  create  an  efficient  and  powerful  computational  algorithm 
over  an  affordable  computer  memory  space.  We  have  considerable  experience  in  the  use 
of  the  conjugate  gradient  algorithm  and  the  FFT  in  both  the  frequency-and  time-domain 
in  our  previous  SBIR  projects.  This  efficiency,  which  does  not  seem  to  be  present  in  any 
other  computational  algorithms,  is  the  basis  of  the  use  of  our  new  higher-order  splines  and 
reciprocal  basis  functions  which  indeed  produces  a  Toeplitz-Hankel  matrix. 
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VVe  will  use  the  conjugate  gradient  together  with  the  fast  Fourier  transform  during 
Phase  II.  We  will  also  address  the  issue  of  the  pre-conditioning  conjugate  gradient.  From  our 
experience,  the  pre-conditioning,  if  used  properly,  can  extensively  accelerate  the  convergence 
rate  of  the  conjugate  gradient  algorithm. 

g)  One  of  the  interesting  features  of  our  newly  developed  efficient  asymptotic  closed  form 
dyadic  Green’s  function  is  that  it  is  amenable  to  transient  electromagnetic  problems  for 
double- layered  anisotropic  slab  (or  half  space).  This  implies  that  an  inverse  fast  Fourier 
transform  of  the  weighted  average  of  our  closed  form  approximation  of  the  Green’s  dyadic 
over  a  proper  frequency  band  can  create  a  correct  response  for  a  prescribed  time  dependent 
electric  point  current  source.  It  should  be  pointed  out,  however,  that  for  the  low  frequency 
components,  where  the  lateral  separation  of  the  source  and  observation  points  are  considered 
to  be  (electrically)  small,  an  efficient  numerical  algorithm  based  on  singularity  removel 
technique,  should  be  used  to  compute  the  Green’s  dyadic. 

We  will  set  up  a  transient  scattering  problem  of  a  perfectly  conducting  strip  within  the 
anisotropic  (uniaxial)  slab  when  a  time  dependent  electric  source  is  located  over  the  slab. 

It  is  noteworthy  that  our  transient  analysis  will  directly  be  applied  for  the  bandwidth 
analysis  of  the  monolithic  microwave  and  millimeter  wave  integrated  circuits  and  the  three 
dimensional  optical  transmission  channels. 

h)  Since  the  geometry  of  the  proposed  model  can  be  very  complex,  we  will  explore  the 
employing  of  a  user-friendly  graphical  pre-  and  post-processor  interface  under  X- windows. 
This  step  is  a  crucial  step  toward  the  commercialization  of  our  microwave  and  millimeter 
wave  integrated  circuit  (MIMIC)  software  product. 

i)  We  have  gained  a  good  deal  of  knowledge  about  the  transputers,  and  other  accelerator 
boards  that  can  be  installed  on  the  386/486  workstation.  These  boards  offer  large-scale 
computing  in  a  workstation  environment.  We  are  also  investigating  the  applicability  of  the 
Intel  Hypercubes  (MIMD),  MasPar  (SIMD),  and  WAVETRACER  (SIMD)  architectures  to 
the  solution  of  our  volume-integral  equation.  The  crucial  feature  of  any  architecture  that 
s  applied  to  our  problem  is  that  it  be  able  to  execute  the  fast  Fourier  transform  (FFT) 
very  rapidly.  We  will  thoroughly  investigate  the  guidelines  that  should  be  followed  in  order 
that  our  computer  aided  engineering  (CAE)  code  will  be  executable  on  massively  parallel 
machines,  as  well  as  on  workstations  with  the  accelerator  boards. 
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-  fi.tlf  +  «,,zzz  ;  |i,  =  fr'tlt  +  /x,iZzz 


Figure  6.1:  Various  bends,  junctions,  discontinuities,  and  couplers  in 
single-layered  anisotropic  uniaxial  slab  used  in  monolithic  microwave 
and  millimeter  wave  integrated  circuit  (MMMIC)  designs 
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fi  ~  +  f«,zzz  ;  hi  =  +  m,ziz 


Figure  6.2:  Finite  micrstrip  phased  array  in  single  as  well  as  double  layered 
grounded  anisotropic  uniaxial  slab 
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*i  —  fi.flt  +  f,,iZ2  ;  Jt,  =  +  H,,zZZ 


Figure  6.3:  Different  bends,  junctions,  and  couplers  of  dielectric  waveg¬ 
uides  in  double-layered  anisotropic  uniaxial  slab  used  in  millimeter  wave 
and  optical  circuit  designs 


77 


Chapter  7 


Appendices 


7.1  Appendix  A:  Eigenvalues  and  Eigenmodes  for  General  Uniaxial 
Media 

In  this  appendix  the  explicit  expression  for  the  eigenmode  na  and  the  eigen  vector-wave 
na{ r,  kt)  for  a  general  uniaxial  anisotropic  material  where  (I  =  etIt  -f  e*zz  ;  p,  = 
Hth  +  MjZz)  will  be  derived. 

> 

As  stated  earlier  in  (2.17),  the  eigenvectors  9a<(p)  corresponding  to  ±Ka  can  be  ex¬ 
pressed  as, 


*S{p)  = 


(  E*’a  *  E*>°’  \  -  (  Zoea  ±  e?i0,  \ 

\  ±HfiC,  +  Hz,o  J  ±h„  +  h*>a  J 


(A.l) 


where  ea,  ha  and  eJ  O,  h,  a  are  normalized  mode- vectors  transverse  and  parallel  to  z  respec¬ 
tively,  for  the  eigenvalues  of  ±k0,  and  Za  is  the  mode  impedance;  and  their  explicit  values 
can  be  found  from  the  source-free  solution  of  (2.7).  Thus,  after  substituting  Vt  — ►  _jk(, 
jL  — ►  -jna  in  (2.14),  I tt  =  €|I4,  and  ptt  =  ptlt  one  will  have, 


wco  Ul<  -  ~~2  —  k(kt|  •  (z  x  ea)  =  Ka  Yaha 

L  u>*(0nofix  j 

uf*o  Pth  -  ~2  -  7  ktkt]  •  (»  x  ha)  =  -Ka  Zaea  , 


(A. 2) 


where  Ya  ~ 

Because  one  is  ultimately  interested  in  the  guided  (^-propagating)  representation  of  the 
vector-wave,  one  can  introduce  a  right  handed  coordinate  basis  of  the  form  (z,  k{,  k"),  where 
the  unit  vector  k{  is  defined  by  the  decomposition  of  k  =  k \kt  +  ZKa  and  k'/  =  z  x  k'(.  Let 
us  now  substitute  an  ansatz  (ej,  hi)  for  q  =  1: 


ei 


K  , 


(A. 3) 
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into  (A.2),  then  one  will  have 


hi  =  z  x  ei  =  k"  , 


u(0et  =  kiYi  , 


-ljuo  ( y.t 


1 

U2€QflQ€z 


k ?) 


«i  Z\  , 


(A.4) 


which  yields  _ _ 

"■ =  ;  z'  =  Ti=  w#,  ’  (A  5) 

where  Atq  =  w2eoMo-  Using  (2.13),  the  normalized  longitudinal  modevector  ez,hz  can  be 
expressed  as, 

e*,i  = - —  ;  /iz.i  =  o  .  (A. 6) 

U€0f, 

As  is  evident  from  (A.6)  the  a  =  1  mode  is  transverse  magnetic  TM  (or  JEJ-mode)  with 
respect  to  the  ^-coordinate. 


For  a  =  2,  one  may  substitute 


into  (A.2), 


which  yields 


e2 

-  k" 

h2 

=  z  x  e2  =  —  k( 

» 

(ft  - 

2  1  *?)  - 

k2  Y2 

fix 

u>H0Ht  = 

«2  Z2 

(A. 7) 


(A.8) 


k2  =  Jklutft  -  —k2t  ;  y2  =  4-  = 
V  ^2 


«2 


upoUt 

Using  (2.13),  the  normalized  longitudinal  modevector  ez,hz  can  be  expressed  as, 

kt 


e*,2  =  0  ;  hZi 2  =  - 


u>hqhz 


(A.9) 


(A- 10) 


As  is  evident  from  (A. 10),  the  a  =  2  mode  is  transverse  electric  TE  (or  iT-mode)  with 
respect  to  the  z-coordinate. 


The  explicit  expressions  for  vector- wave  modes  of  a  general  uniaxail  anisotropic  medium 
can  now  be  obtained  (see  (2.17)  after  incorporating  (A.3),  (A.5),  (A.6),  (A. 7),  (A.9),  and 
(A. 10)  into  (A.l). 
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7.2  Appendix  B:  Reciprocity  and  Orthogonality  Relations  in  Anisotropic 
Media 

In  this  appendix  a  set  of  reciprocity  and  orthogonality  relations  for  anisotropic  media  will 
be  discussed.  These  relations  are  useful  in  the  derivation  of  the  dyadic  Green’s  function  for 
unbounded  as  well  as  multi-layered  uniaxial  anisotropic  media. 

The  Maxwell  curl  equations  for  a  gyroelectric  source  excited  medium  can  be  expressed 
as. 


V  x  E  =  - jw/io/xH  -  M  , 

VxH  =  jwe0I  ■  E  +  J  ,  (B.l) 

where  (^oA>fo<0  arc  the  constitutive  tensors  of  the  medium,  and  J  ?:id  M  are  electric  and 
magnetic  current  sources,  respectively.  One  can  introduce  an  associated  adjoint  medium 
with  the  constitutive  parameters  of  (^oA  +  ifo^+),  where  in  the  adjoint  medium  the  Maxwell 
curl  equations  can  be  expressed  as, 

V  x  E+  =  -jwji0A+H+  -  M+ 

V  x  H+  -  ju e0e+  •  E+  +  J+  .  (B.2) 

The  appropriate  reciprocity  relation  can  be  sought  by  calculating  V  •  (E  x  H+)  = 
H+  •  V  x  E  -  E  ■  V  x  H+  and  V  •  (E+  X  H),  from  (B.l)  and  (B.2)  taking  their  difference: 

V  •  (E  x  H+  -  E+  x  H)  =  +  jueo  (e+  •  I  •  E  —  E  •  f+  •  Ef ) 

(H+  •  A  H  -  H  A+  H+)  (B.3) 

H+  •  M  +  E+  •  J  +  H  •  M+  -  E  •  J+  . 

If  the  transpose  of  the  permiability  and  permittivity  tensors  of  the  medium  are  the  same 
as  the  permiability  and  permittivity  tensors  of  the  adjoint  mtuium,  (A  =  A+  ;  I  =  e+, 
where  (~)  denotes  transpose),  then  one  can  identify, 

E+^T-E  =  E  • !  •  E+  -  E-S+-E-*  , 

H+^A  H  =  H  A  H+  =  H  A+  H+  , 

and  the  first  four  terms  of  the  right  hand  side  will  cancel.  Upon  introducing  a  closed  volume 
V  bounded  by  the  surface  5,  one  will  have 

ds  ■  (E  x  H+  -  E+  x  H)  =  J  dv  (-H+  •  M  +  E+  •  J  +  H  •  M+  -  E-  J+)  .  (B.4) 

Equation  (B.4X  is  the  r  ciprocity  relation  in  an  anisotropic  (symmetric  gyroelectric)  medium 
wherein  {ft  =  ji  ;  e  =  e),  and  tht-  adjoint  medium  will  be  the  same  as  its  associated  medium. 
A  general  uniaxial  medium  is  an  special  case  of  symmetric  gyroelectric  medium. 
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The  orthogonality  of  the  z-guided  vector-wave  modes 


n<(k„ r)  =  £ (  ih'j+h’;:  ) <B5> 

can  be  inferred  by  substituting  for  the  divergence  operation,  V-  =  (vt  +  •  ,  in  (B.3) 

and  applying  the  surface  divergence  theorem  to  a  pair  of  source  free  fields  of  /7„(kt,  r)  and 

< 

Z7^(-k{,r)  [22,  pages  174-179];  this  yields, 

J  ds  Vt  •  (±Zaea( kt)  x  h^k*)  -  Z0e0(k[)  x  h«(kt)) 

=  j>  Mb-  (±Zaea{kt)  x  Mk{)  -  Z0e0(K)  x  hQ(kt))  (B.6) 

=  YzJs  dsi  ‘  (±Z“e“(k‘)  x  Mk'*)  -  Wk;)  X  h0(kt))  . 

The  contour  integral  vanishes  because  of  the  radiation  condition  as  p  — ►  oo1.  Thus  one  will 
have 

dsi  •  (±Zaea(kt)  x  h,i{k't)  -  Z0e0(k')  x  ha(k())  e--*’<k'-k!>  =  0  ,  (B.7) 
Js 

which  yields, 

hj  ■  (z  x  ea)  -  ^a/3  ,  (B.8) 


®/3  '  (fia  X  z)  —  ^ar/3 


(B.9) 


It  is  noteworthy  that  the  orthogonality  relations  stated  here  do  not  distinguish  between 
the  modes  for  +«„  and  -Ka,  which  is,  of  course,  the  direct  result  of  the  reflection  symmetry 
property  of  the  medium. 


‘The  surface  divergence  theorem  is  applied  to  the  open  surface  5  here  with  the  radiation  condition 
implied  as  p  — *  oo. 


81 


7.3 


Appendix  C:  Essential  expressions  related  to  g£{£,  z,  z')  and  g-fit,  z,  z ') 


,> 

In  this  appendix  explicit  expressions  for  gif(£,  z,z')  and  g-<(£,  z,  z')  along  with  the  the 

residues  of  the  surface  wave  and  leaky  wave  poles,  which  are  related  to  the  roots  of  De 

> 

and  Du  are  given.  Also,  expressions  for  iV^(l),  lim£_i 

H 


H 


,  £>£(1),  and 
H 


lim^_i  —  £2  which  are  employed  in  Chapter  4  to  construct  the  uniform 

asymptotic  representation  of  the  dyadic  Green’s  function  are  presented  in  this  appendix. 


We  first  start  with  TM  or  (E-mode)  terms,  gjf(£,  z,  z'),  and  subsequently  the  TE  or  (H- 

„>  >  > 

mode)  terms,  g^it,  z,  z'),  will  be  presented.  The  functions  <?'<(£,  z,  z'),  and  z,  z'), 

can  be  expressed  as 


De(  0 

N>jt„U,z,z') 

Dh(0 


(C.1) 

(C.2) 


where  functions  NjjtE,  De,  Ntt,H ,  and  Djf  will  be  defined  shortly.  Between  gf},  and  gf}  the 
following  relations  hold: 


g\f{z,z')  = 

9u{z',z)  , 

(C.3) 

9^,z')  = 

(C.4) 

w 

** 

II 

1 

*  A 

(C.5) 

9u  (*>*')  = 

** 

Ah 
^  ■** 

1 

(C.6) 

■-97{z',z)  . 

(C.7) 

The  functions  N^e,  and  De  (in  view  of  (2.18),  (2.36),  and  (2.49))  are  defined  as 


N>e(z,z')  = 


e— r«2(*— *')(!  _ 

1  +  e-]2K& 


[Ui  (1  +  e-j2'^irrr))  +  U2(  1  -  e-^t*1-*’)] 

(C.8) 
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N>e(z,  z')  =  [t;,  (1  +  e-ia-4(ST=D)  +  (1  -  e-i2«i(IT=7))] 

1  +  e  J  2  '  (C.9) 


#tU(*’Z,)  =  - 


1  +  e-j2«' Js 


[tfi  (1  -  e-22"i(*i-z'))  +  l/2(H-  e-ia"4(*i-*'))l 

(C.10) 


r>  ./%  _  e-2'ti(z-*')(l  +  e-J2',i(*'-**)) 


K,e(*,z')  = 


1  4.  e-j2«' J2 


[tfj  (1  -  e-j2*i(^I))  +u2(l  +  e-J2«i(^D)] 

(C.ll) 


and 


where, 


s*(0  =  ^(0 +  i  tan KjrfaCWO  , 


17i(0  =  ^1  {Zo  ■+•  j  Z%  tan  K[di) 


U2(£)  =  Z2  (Z\  +  j  Zq  tan  Kjdi)  , 
and  parameters  z,  —  Zj,  d,,  /c',  are  defined  as 


(C.12) 


(C.13) 


-  Zj  =  ko(zi-Zj)  ;  d,  =  k0di  ;  k-  =  Jn,  -  n'f  2  ;  n<  =  m,t€i,t  ;  K  =  —  ,  (C.14) 

e»-* 

and  the  normalized  mode  impedance  Z,  is  defined  as 


Z.  =  ^  . 
e.,t 


(C.15) 


Note  that  the  parameters  and  Z,  are  respectively  normalized  to  ho  and  u;eo,  and  thus 
n0  =  n'0  =  1. 

,> 

The  residue  of  g£(z,  z ')  at  £  =  £p>  which  is  associated  with  the  TM  (or  E-mode)  surface 
wave  or  leaky  wave  mode  is  given  by 


Res  g 


= 

'  U=<„/ 


(c.i6) 
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where 


JlD *<« 


■A 


Zi 


1  ,n\  ( tan  K\d\  d 

-  +  J-  — r  + 


«0 


COS2(Kjdi) 


+  — ^7"  (^o  +  j  Z\  tan  Kj d\) 

fiKi 

i  f  tan  K[di  n'jKo  d\  \ 

\  Kq  k[  COS2  KjJi  J 

(tan  4d,  +  _*_)  (Zl  +  iZo  tan  sWl)l )  . 
COS2(K'2a2)  /  J  J 


+j 


Z2  tan  «2</2 


ni 


(C.17) 


+  -^ 
^2 


Also,  explicit  expressions  for  N^E(1),  lim^_x  [\/l  —  Is  i^E(0],DE( l),andlim^x  [v/T^T 
are  given  by 


N>e(1,z,z')  = 


-i\Jn2  ~  _  e-j2yjn2  -  n'2(z‘- 


*2) 


) 


1  +  e 


-j2  yjn2  —  n  2  <f2 


tan  y/nx  -  n'idxj  (1  +  e~lyfi2  (*1  z))  + 

(1_  e~i2  \Jn2  ~  n’2(—)^ 


yjn2  —  n'2  ^ 

/nx  -  ni 

€ 2 

fl 

(C.18) 


N>e(1,z,z')  = 


-j\Jn2  -  n'2  (*~*')(1  _  e->2 yjn2  -  n'2(z'-z2) 


) 


.nx  —  ni 

lj— tan 


1  e-j2y/n2  -  n'2i- 


-  n'2  \Jn\-  n'j 

e2  fi 


(i+e-Vn2-n2(irr7))  + 

^  _  e~22\Jn2  -  rc2(*l-*)^ 


(C.19) 


iV<£(l,2,/)=  - 


0-i\Jn2  -  n'2(z-z')^  _  e-j2yjn2  -  n'2(*'- 


1  +  e 


-j2yjn2  -  n'2d2 


(C.20) 


i^(o] 
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(C.27) 


Also  the  following  relations  exists 


(C.28) 


lim 


N^c(l,z',z)  , 


(C.29) 
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N>,B(l,z,z')  =  ~N<'E(1,z',z) 
lim  v/l-€2  !*>£(*, 2')  =  ~  g*  JcNuAz'' z)  > 


(C.30) 


N>E(l,z,z')  =  -AT<£(l,2',z) 

Um  v/i  -  £2  4 XZ'E(z,z')  =  -  lim  ^/l  -£2  ^Nft  E{z',  z) 


(C.31) 


As  mentioned  earlier,  the  T.E  or  (H-mode)  terms,  g"f(£,  z,z'),  is  given  by 


ff”/  (£,  A  = 


,x  N?<h(£,z,  z q 

j  0»(O  ’ 


(C.32) 


where  the  functions  and  DE  (in  view  of  (2.20),  (2.36),  and  (2.49))  are  defined  as 


JV>„  =  [v,  (I  -  ,-»»s"(3r=))  +  F,  (1  +  «-»«!(!r=^)]  , 


where 


•0//(0  =  *2(f)  +i  tanK2d2  Vx(£)  , 


Vi(0  =  Vi  (Yo  +  jYt  tan^dx)  , 


Va«)  =  Y2(Y1+jY0t^K’{dl)  , 
and  the  parameter  k "  is  defined  as 

<  =  ;  n"  =  —  , 

M*,* 

and  the  normalized  mode  admittance  y;  is  given  by 


(C.33) 

(C.34) 


(C.35) 


(C.36) 


y,  =  -2-  . 

mm 


(C.37) 
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Note  that  the  parameters  /c,-,  and  Yi  are  respectively  normalized  to  ko,  and  and  thus 
n0  =  ti'q  =  1. 

> 

The  residue  of  {/"fiz,  z')  at  £  =  <fp»  which  is  associated  with  the  T  E  (or  H-mode)  surface 
wave  or  leaky  wave  mode  is  given  by 


Res  2') 


'-v  "  ' 


where 

£ 

d£ 


[  [Mi*i  \  «o  «"  cos2(«'1d1; 


H - ^77  (Yi  +  j  Yo  tan  n"di) 

H2K2 


+j 

JjZ 


,,  1  ,  .n'/  / tan  /ci'di  di  \ 

yt  tan  —  +  J—  - — ~  + - 

[«o  \  Ki  cos2  KWi  / 

-  tan  n2d2  + 


MlKl 


V  cosV2'd2))  {Yo  +  jYi  tan 


(C.38) 


(C.39) 


*»>]}  • 


Also,  explicit  expressions  for  Nfi H(  1),  lim^_i  —  !*  Dh(1),  ^d  limf_,  [^/W^wtO] 

are  given  by 


N tt.H (^’  — 


~i\J n2  ~  _  e-j2y/n2  -  n%(z'-z2)] 


1  +  e 


-j2  yjn2  ~  n%J7 


_  _//  , 

.Tl\  Hi  / 

J-7-— —  tan  V"1  “  "tdl 


(Ml)2 


_  — 2 \J n2  ~  n'2 


(1  -  e 


)  + 


(C-40) 


y/»2  ~  7*2  y^i  -  n” 
M2  Mi 


(1  +  e 


-ji\Jn2  -  n'2'(7T=7) 


lim 

<-1 


->\,/n2  -  "2  (*"*')( i  _  e~j2\in 2  '  n2,(*'-**)1 


1  +  e 


-j2 


yj^2  —  n2  dj 


'J™1  _  e~J\Jn2  -  n2  (Zi-z) 


Mi 


)  + 
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7.4  Appendix  D:  Limiting  values  of  g^i^z.z')  and  g”f(£,z,z’)  as  £  — *  oo 

<>  „> 

In  this  appendix  the  limmiting  values  related  to  gif(z,z')  and  gi^(z,zl)  as  £  — *  oo  which 
are  being  used  in  (3.33)  will  be  given. 

Note  that 


lim  k' 

(~>ao 

lim  k" 
€—<» 


Hence  (see  (C.l)  -  (C.ll)  and  (C.32)  -  (C.39)) 


(D.l) 


lim  gt?(z,z') 
*oo 


_ _ 2/K  _  ^5  _ _ 

(1  ~  a 


_ _ _ _ 

lim  g'^(z,z')  ->  e-«y^i(*-*')  -  ~^= - e~(^2xi -(*+*'»  , 

<»  y"'i  ,  y  n2 

e,  +  e2 

lim  g  'tfiz,  z')  — ►  e-fV^(z-z^  +  — - - -  e~f\/"?(2*l~(z+i')) 

#*1  ~  « 


Also,  as  £  — »  00  the  following  relations  exit: 


(D.2) 


(D-3) 


(D.4) 


gtf(z,zf) 

=  ~9u  {*',*)  , 

(D.5) 

=  9tf(z,z')  , 

(D.6) 

9,f(z,z') 

=  -9tf(z\z)  , 

(D.7) 

9'tf(z,z') 

=  -$£(*',*)  , 

(D.8) 
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and 


&*(*>*')  = 

9*(z,z>)  , 

(D.9) 

<C>» 

Ra 

II 

feV>z)  , 

(D.10) 

$”*(*>  z')  - 

=  $«>(*,  *0  • 

(D.ll) 
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